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Preface to these three texts 

The three texts in this one cover, entitled ‘First-order partial differential equations’ (Part I), ‘Partial differential equations: 
classification and canonical forms’ (Part II) and ‘Partial differential equations: method of separation of variables and 
similarity & travelling-wave solutions’ (Part III), are three of the ‘Notebook’ series available as additional and background 
reading to students at Newcastle University (UK). These three together present an introduction to all the ideas that are 
usually met in a fairly comprehensive study of partial differential equations, as encountered by applied mathematicians at 
university level. The material in some of Part I, and also some of Part II, is likely to be that encountered by all students; 
the rest of the material expands on this, going both further and deeper. The aim, therefore, has been to present the 
standard ideas on a broader canvas (but as relevant to the methods employed in applied mathematics), and to show how 
the subject can be developed. All the familiar topics are here, but the text is intended, primarily, to broaden and expand 
the experience of those who already have some knowledge and interest in the subject. 

Each text is designed to be equivalent to a traditional text, or part of a text, which covers the relevant material, but in a way 
that moves beyond an elementary discussion. The development is based on careful derivations and descriptions, supported 
by many worked examples and a few set exercises (with answers provided). The necessary background is described in the 
preface to each Part, and there is a comprehensive index, covering the three parts, at the end. 
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List of examples 

Here is a list of the various examples that are discussed in this text. 

'y 

General solution of 2yu x + uu y = 2yu .p.16 

Show u(x,y ) = e- v F[(l + xw)e - v ] is a solution of 2 yu x +uu y = 2 yu .p.16 

General solution of xu x + 2 xuu y = u and solutions with: (a) u -2x on y = 2x + 1; 

2 3 2 3 

(b) u - 2x on y = 3x ; (c) u = x on y = x -1.p.20 

General solution of u x +2xu y = u .p.22 

General solution of yu x - xu y = 2 xyu .p.23 

General solution of u x + e x u y + e z u z - (2x + e' T )e" =0.p.23 

Solve ODEs that describe u x u y = u i.e. pq-u .p.32 

Solution of u x u y = u with u-t on x-t, y = 1 +1 .p.34 

Solution of u x + u y = 1 with u - At on x = y-t .p.35 

Find complete integral of u x u y =u .p.41 

9 9 

Find complete integral, general solution, singular solution of u x +u y =1 + 2 u .p.42 

Solution of xu x n y + yu y = 1 with u-t on x = 2 1, y = 0.p.44 

9 9 

Find complete integral of y(u x -u y ) + uu y = 0 and solutions with 

2 

(a) u = 3t on x = 2t, y = t;( b) u - 2t on x-t , y = 0.p.45 
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Preface 


Preface 


This text is intended to provide an introduction to the standard methods that are used for the solution of first-order partial 
differential equations. Some of these ideas are likely to be introduced, probably in a course on mathematical methods 
during the second year of a degree programme with, perhaps, more detail in a third year. The material has been written 
to provide a general - but broad - introduction to the relevant ideas, and not as a text closely linked to a specific module 
or course of study. Indeed, the intention is to present the material so that it can be used as an adjunct to a number of 
different courses - or simply to help the reader gain a deeper understanding of these important techniques. The aim is 
to go beyond the methods and details that are presented in a conventional study, but all the standard ideas are discussed 
here (and can be accessed through the comprehensive index). 

It is assumed that the reader has a basic knowledge of, and practical experience in, the methods for solving elementary 
ordinary differential equations, typically studied in the first year of a mathematics (or physics or engineering) programme. 
However, the development of the relevant and important geometrical ideas is not assumed; these are carefully described 
here. This brief text does not attempt to include any detailed, ‘physical’ applications of these equations; this is properly left 
to a specific module that might be offered in a conventional applied mathematics or engineering or physics programme. 
However, a few important examples of these equations will be included, which relate to specific areas of (applied) 
mathematical interest. 

The approach that we adopt is to present some general ideas, which might involve a notation, or a definition, or a method 
of solution, but most particularly detailed applications of the ideas explained through a number of carefully worked 
examples - we present 13. A small number of exercises, with answers, are also offered, although it must be emphasised 
that this notebook is not designed to be a comprehensive text in the conventional sense. 
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1 Introduction 


The study of partial differential equations (PDEs), both first and second order, has a long and illustrious history. In the 
very early days, second order equations received the greater attention (essentially because they appeared more naturally 
and directly in problems with a physical basis). It was in the 1770s that Lagrange turned to the problem of solving first 
order PDEs. [J.-L. Lagrange, 1736-1813, an Italian-born, French mathematician, made contributions to the theory of 
functions, to the solution of various types of equations, to number theory and to the calculus of variations. However, his 
most significant work was on analytical mechanics.] 


All the essential results (in two independent variables) were developed by Lagrange, although Clairaut (in 1739) and 
dAlembert (in 1744) had considered some simpler, special PDEs. The bulk of what we describe here is based on Lagrange’s 
ideas, although important interpretations and generalisations were added by Monge, in particular. The extension to 
higher dimensions was completed by Cauchy (in 1819), and his ‘method of characteristics’ is often the terminology used 
to describe Lagrange’s approach. 


1.1 Types of equation 


We shall concern ourselves with single equations in one unknown - so we exclude, for example, two coupled equations 
in two unknowns - and, almost exclusively, in two independent variables. Thus we shall seek solutions, u{x,y ) , of first 
order PDEs, expressed most generally as 


du dll 

f(x,y,u,—,—) = 0, 
ox dy 


although we shall mention, briefly, the corresponding problem for u{X\,X 2,....,X„), where appropriate. It will be 
convenient, however, to start with a simpler, general equation: 


a{x,y,u ) 


du 

dx 


+ b(x,y,u) 


du 

dy 


c(x,y,u), 


the quasi-linear equation i.e. one that is linear in the two first partial derivatives. This case includes the linear first order PDE: 

a(x,y)^- + b(x,y)^- = c(x,y)u , 
dx dy 


and the semi-linear equation 


a(x, y) ^- + b(x,y)^- = c(x, y,u), 
dx dy 


each for suitable, given functions a, b and c. 
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The quasi-linear structure - linear in ll x and Uy , using the shorthand notation - leads to a fairly straightforward method 
of solution. This is generalised and extended in order to solve the most general PDE of this type, which involves U x and/ 
or Uy no longer of degree 1. Thus we shall discuss the solution of equations that may look like 

_ du du 

2 y -1- u — 

dx dy 



or 


X 


fdu} 

\dxj 


2 k a A 2 


+ 


du 


= 2 xyu . 


As we shall see, the description and construction of the solution requires some relevant geometrical ideas. 


Exercises 1 

Classify these equations as linear, semi-linear, quasi-linear or general: 

7 2 3 2 

(a) ( v-u)u x + xUy = xv + u ; (b) uu~ -xu y — —u ; (c) x u x +(y-x)u,, = vsinw; 
(d) (sin y)u x -c x u y = Q y u ; (e) u x +sin(w > ,) = u . 


sfsfX-X-X-X-X-X-X-X-X-X-X-X-X-X- 

x-x-x-x-x-x-x-x-x-x- 
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2 The quasi-linear equation 

We develop first, with care, the solution of the most general equation written in the form 

a(x,y,u)u x +b{x,y,u)u y =c(x,y,u), 

where we have used subscripts to denote the partial derivatives. The coefficients a, b and c need not be analytic functions 
of their arguments for the most general considerations. However here, in order to present the conventional and complete 
theory, we shall assume that each of a, b and c is in the class of C ' functions i.e. those that possess continuous first partial 
derivatives in all three arguments (at least, in some domain D). Further, we will normally aim to seek solutions subject to 
a given boundary condition, namely, that u{x, y) is prescribed on a known curve in the ( x , y)-plane (which will need 
to be in D for a solution to exist); this is usually called the Cauchy problem. 

2.1 Of surfaces and tangents 

Let us suppose that we have a solution U =li(x,y ) which is represented by the surface z = u(x,y) in Cartesian 3-space. 
Now it is a familiar (and readily derived) property that the surface u(x, y) — Z = 0 has a normal vector (it x ,U y ,— 1) 
at every point on the surface. Further, we introduce the vector ( a,b,C ), where a, b and c are the given coefficients of 
the quasi-linear equation i.e. 


au x +bu y -c = 0 becomes (a,b,c)-(u x ,u y ,-\) = 0 . 

Thus the vectors ( a,b,C ) and (u x , U y , — 1) are orthogonal, and because (u x ,U y ,— 1) is normal to the surface 
ll(x, y) — Z = 0 , so Z = u(x, y) is in the tangent plane defined by the vector (a, b, c). In other words, the surface 
which is a solution, z = w(x, V) , is the surface which has ( a,b,C ) as its tangent everywhere: the tangent plane. 

The solution of the quasi-linear equation can therefore be expressed by the description of the tangent plane in terms of 
the slope of this surface: 

d z _ c dz _ c 

d x a d y b 

However, the required solution (surface) is Z = ll(x, y) , and so we may write 

du c(x,v,u) d u c(x,y,u) 

-=- 1 - and -=---, 

dx a(x,y,u ) dy b(x,y,u ) 

or, equivalently, one of this pair together with 

dv _ b(x,y,u ) 
dx a(x,y,u ) 
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This last equation defines a family of curves (but dependent on u ) that sit in the solution surface, and are usually called 
characteristics (after Cauchy); the set of equations is usually called the characteristic equations of the PDE. The result of the 
integration of these two (coupled) ordinary differential equations (ODEs) is a two-parameter family - the two arbitrary 
constants of integration - and then a general solution can be obtained by invoking a general functional relation between 
the two constants. This is sufficient information, at this stage, to enable us to discuss an example. 

Example 1 

2 

Find a general solution of 2yu x + till., = 2yiC. 


d u 2 yu 

First we have the ODE -= —-= 1 Vli (for U A (J ), which gives 

dy U 

J—= 2j vdv and so ln|w| = y 2 + constant i.e. u = Ae y , 

VI 

2 

dy u Ae y 

where A is an arbitrary constant. Now consider — =-=-, then 

dx 2 v 2 v 

2 J ve _v dv = A J dx or -e~ y = Ax-B i.e. Ax + e~ y = B , 

where B is the second arbitrary constant. The general solution is expressed by writing A = F(B ) , where 


B = Ax + e y = (1 + x«)e 


which gives 


u(x,y) = e } i 7 [(l + xn)e ], 


an implicit relation for u{x,y ) , where F(.) is an arbitrary function. Note that the zero solution is recovered by setting 

F = 0. 


Comment: That we have obtained an implicit, rather than explicit, representation of the solution is to be expected: the 
original PDE is nonlinear. In the light of this complication, it is a useful exercise to confirm, by direct substitution, that 
we do indeed have a solution of the equation for arbitrary F; this requires a little care. 

Example 2 

2 _ 2 2 

Show that U (x, v) = e v /*’[( 1 + XII )c y ] is a solution of 2yu x F ully = 2yu , for arbitrary F. 
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First, we see that U x = e - ' ( Xll x +u)e •’ F'(^) (where £ = (1 + xw)e •' ), which gives 

uF' 


u x = 


1 - xF' 


(for 1 — xF' ^0); 


similarly, we obtain U v = 2 ye- F + C' \xu y Q - — 2 V (1 + XU )c ' ]F , which gives 


11 y = 


~y _2 i.n -l v„WT 177' 

2ve- v F-2 v(1 + xt7)F' 


1-xF' 


Now we form 


2v»x + uu y 


2vuF' + 2yue y F-2yu(\ + xu)F' 
1 -xF' 

2 yue y F -2xyu 2 F' 

1 -xF' 


2yuQ y F -2xyue y FF' 
1 -xF' 


2yue y2 F(\-xF') 
1 -xF' 


2yue y F = 2yu 2 , 



AACSB 

. ACCREDITED , 


CLICK HERE 

to discover why both socially 
and academically the University 
of Groningen is one ofthe best 
places for a student to be 


www.rug.nl/feb/education 


Excellent Economics and Business programmes at: 

university of 
groningen 


“The perfect start 
of a successful, 
international career.” 


17 

Download free eBooks at bookboon.com 











An introduction to partial differential equations 


The quasi-linear equation 


which recovers the original equation: we do indeed have a general solution (for arbitrary, differentiable F, and where 

1-jcFVO). 


In this example, we observe that the characteristic lines are 

_ 2 

(1 + xw)e y = constant, 

2 

along which we •’ is a constant. We now examine more carefully the connection between the characteristic lines and 
the curve on which specific initial data is prescribed. 

2.2 The Cauchy (or initial value) problem 

The most convenient way to proceed is to introduce a parametric representation of the various curves in the problem - 
both those required for the general solution and that associated with the initial data. First, let the initial data (i.e. u given 
on a prescribed curve) be expressed as 


u=u(t ) on x = x(t), y = y(t), 

where t is a parameter that maps out the curve, and the u on it. Of course, this is equivalent to stating that II (x, y ) 
is given on a curve, g(x, v) = 0, say. Now we turn our attention to the characteristic lines, and the general solution 
defined on them. 

The pair of ODEs that describe the solution can be recast in the symmetric form 

dx _ dy _ du 
a(x, y,u) b(x,y,u ) c(x,y,u)' 

and then any two pairings produce an appropriate pair of ODEs (as in §2.1). However, we may extend this variant of the 
equations by introducing the parameter, s, defined by the construction 

dx _ dy _ dw _ ^ 
a b c 


in other words, we may write 

dx . dv , , dw . . 

— = a(x, y, u),— = b(x, y,u), — = c(x, y, u). 
ds as ds 

These three (coupled) ODEs are the parametric equivalent to the pair of equations derived in §2.1. A solution of this set 
will then be expressed as 


w = u{s) on x = x(s), v = t(s) , 
which mirrors the description of the initial data. 
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Now a solution that satisfies the PDE and the initial data will necessarily depend on the two parameters: t chooses a point 
on the initial-data curve, and then s moves the solution away from this curve along a characteristic; see the sketch below. 



Thus a complete solution, which satisfies given initial data, will depend on the two parameters, s and t: 

u = u(s,t ) on x = x(s,t), y = v(s,t). 

The only possible complication that might arise is when we attempt to construct u(x,y) (which, presumably, is what 
we are seeking). 


To accomplish this construction, in principle, we must solve X = x(s,t) and 


y = y(s,t ) to find S = s(x, y) and t = t(x, y) ; these are then used to obtain u(x, y) . The solution for s and t exists 
provided that the Jacobian of the transformation is non-zero i.e. 


j_ d ( x ,y) 

d(s,t) 


x f 


y s 

y t 


= x s y t - x ty s * °- 


If J = 0 anywhere on the initial-data curve, then the solution fails (for this particular data): the solution does not exist. 
So, provided J ^ 0, a (unique) solution satisfying the PDE and the initial data exists; but what if we do find that J — 0 ? 


If J = 0 then 
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let us construct 


and so 


Similarly, we can obtain 


that is, if J = 0 , we have 


w = x ty s 


or 


y t 


X 


s 


y s 


a 

V 


du 

d t 


—u(x,y) = u x x t +u y y t 


= u x^Y +u y yt=^ ( au x + bu y )y t 


U t _ au x + bu y y t _ au x +bu y y f _ y f 
c c b au x + bu y b b 


u t _ au x +bu y X{ _ au x +bu y Xf _ x t 
c c a au x + bu y a a 


^ L = ^ L = y L 

cab 


du _dx _ d y 

which is equivalent to the original equations i.e. -— — — ——. Thus the initial data, in this case, must itself be a 

cab 

solution of the original PDE; in other words, the curve on which u is prescribed is one of the characteristics, and the u 
on it must satisfy the PDE. The solution is therefore completely determined on this particular curve (characteristic line), 
but is otherwise non-unique on all the other characteristics. These various points are included in the next example. 


Example 3 

Find the general solution of XU x + 2xuu y = U and then seek those solutions (if they exist) subject to: (a) ll = 2x on 
y = 2x 2 + 1 ; (b) U= 2x 2 on y = 3x 3 ; (c) U = X 2 on y = X 3 -1. 


The equation is equivalent to the system 


dx dv du du u dy _ 

— = — 1 — = — e.g. — = — and = 2 u . 

x 2 xu u dx x dx 


fd« rdx 

Thus -= — and so U = Ax (A an arbitrary constant), 

J U J X 

dv 2 

and then = 2it = 2 Ax which gives V = Ax + B (B the second arbitrary constant), 
dx 

2 

The general solution is described by A = F(B) , with B = y — Ax = y — XU , i.e. 

w(x, y) = xF(B ) = xF(y-xu). 


20 


Download free eBooks at bookboon.com 








An introduction to partial differential equations 


The quasi-linear equation 


2 2 2 

(a) This requires 2x = F(2x +1 — 2x ) = xF(l) , so that A(l) = 2 : the solution is U = 2x on y = 2x“ + 1, and 

non-unique elsewhere. 

2 3 3 3 1/3 

(b) Now we have 2x = xF(3x — 2x ) = xF(x ) , and so F(z) = 2 z ' ; the solution is then 

u(x, y ) = 2x( v - xu) 1 / 3 , 

which can be recast as a cubic, and analysed further (if appropriate): 

u 3 +8x 4 m = 8 vx 3 . 
x 2 = xF(x 3 -\-x 3 } F(-l) = x 

(c) Finally, we require v ’ or v ' , and this is clearly impossible (for x is an independent 

variable): no solution exists that satisfies this particular initial data. 


Comment: The characteristic lines here are y — XU = constant which, as we see, depend on u. In (a), with U = 2x 

2 2 

, the characteristic lines become V — 2x“ = constant , and U = 2x with V = 2x + constant is a solution of the 

original equation. In (b), the characteristic lines evaluated for the initial data are V — 2x = constant , but the data is 

3 3 

on y — 3x = 0 . In the last exercise, (c), the characteristic lines for this initial data become y — X = constant , and 

3 2 

the data is given on V — X = 1, but with U = X“ this is not a solution of the original PDE. 
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2.3 The semi-linear and linear equations 


The semi-linear and linear equations are described by 

a(x,y)u x +b(x,y)u y 


f c(x, y, u), semi-linear 
[c(x,y)u, linear 


for which the underlying system of ODEs is, correspondingly, 

dx _ d y fd u/c(c,y,u) 

a(x,y ) b(x,y ) \du/uc(x,y). 

The first ODE, in both cases, is 

dv _ b(x, y) 
dx a(x,v) 

which describes the characteristic lines, now defined independently of the solution u(x, v) . The second equation of the 
pair is, for example, 


du c(x,y,u) . x Id u c(x,y) n wi . 

-=- (semi-linear) or-=- (U U) (linear), 

dx a(x,y) u dx a(x, y) 


and hence we may expect the construction of the solutions in these two cases to be particularly straightforward. 


Example 4 

2 

Find the general solution of U x + 2xUy = U . 


, , dx dv du , , 

Here we have the set — = = —— , which we choose to write as the pair 


1 2x 


dv „ du ? 

= 2x and -= U “ ; 

dx dx 


2 1 1 

these are integrated directly to give y = X + A and-= X — B or U — ——- , respectively, where A and B are the 

arbitrary constants of integration. The general solution is then recovered from B = F(A ) i.e 


u(x,y) = 


1 


F(y-x 2 )-x 


for the arbitrary function F(.). 


2 

In this example, the characteristic lines are y — X = constant; for reference, which we shall recall below, we observe 
the way in which the arbitrary function appears in this solution. Now we consider an example of a linear equation. 
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Example 5 

Find the general solution of VII x — XU y = 2XVLI. 


dx dv du dy xdu 

The system of ODEs is — = — = - e.g. — = - and -= Ixu (y A (1 ); these give, respectively, 

y -x 2 xyu dx v dx 

y 1 + x 2 = A and In |z/| = x 2 + constant or u = Be x , 
where A and B are the arbitrary constants of integration. The general solution is then given by the choice B F(A ) : 


u(x,y) = e x F{x 2 Ay 1 ). 


Comment: This linear PDE possesses a solution in which the arbitrary function in the general solution appears in the 
form U cc F i.e. u is linear in F; this is a necessary property of the linear PDE. This is to be compared with the way in 
which the arbitrary function appears in the solution constructed in Example 4 - this equation was nonlinear. 


2.4 The quasi-linear equation in n independent variables 

The development described in §2.1 is readily extended to n independent variables, X (i = 1 , 2 ,...n ). The quasi-linear 
equation in n independent variables is written as 

du 


Y j a i (x l ,X2,...,x n ,u) — = c(x l ,x 2 ,...,x n ,u), 


i =1 


dx,- 


and then the slopes of the tangents on the «-dimensional solution surface, a surface sitting in n +1 dimensions, generate 
the system of ODEs: 

dxj _ dx 2 _ _ dx„ _ du 

Cl 2 Clyj C 


It is clear that the formulation we gave for the 2-D surface in 3-space goes over to higher dimensions in the obvious (and 
rather neat) way detailed above. The general solution is then described by a general functional relation between the n 
arbitrary constants of integration; we show this in the next example. 


Example 6 

Find the general solution of U x + Q X U y +e~U z — ( 2 x + = 0 . 


23 


Download free eBooks at bookboon.com 









An introduction to partial differential equations 


The quasi-linear equation 


The system of ODEs is 


dr dv dz d u 


1 e x e z (2 x + e x )u 


djf i ^ z , dw i 

e.g. — = e ,— = e and— = (2r + e )e 

dr dr dr 


The solutions are, respectively, 


y-e x = A,x + e " = B, e “ + x 2 + e x = C■, 


the general solution is then obtained by writing C = F(A,B) : 

—u p/ x , —z\ 2 x 

e = F(y-e ,r + e )-r -e , 


from which u follows directly by taking logarithms; the function is arbitrary. 


We have explored the simplest type of first order PDEs: those linear in the first partial derivatives. It should come as no 
surprise that, if the PDE is not linear in both U x and U v , then complications are to be expected. In the next chapter we 
investigate, in detail, the nature and construction of solutions of the most general, first order PDE in two independent 
variables. 

Exercises 2 

1 . Find the general solution of (1 — xu)u x + v(2r + u)u y = 2x( 1 — Xll ) , and then that solution which satisfies 
U = Q y on r = 0 ■ 

2. Find the general solution of e - ' U x + XU y = XU ~ , and then that solution which satisfies It = exp(r“) on V = 0. 

3. Find the general solution of U x —2 XUU y = 0, and then the solutions (if they exist) that satisfy (a) U = X ' on 

-J ^ 

y = 2r; (b) u = x on y = r . 


x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x- 

********** 
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3 The general equation 


This equation is written in the form 


/ 


x,y,u, 


du du 
dx ’ dry 


= 0 , 


which is conventionally expressed in the notation 


f(x,y,u,p,q) = 0 


du 

p=—,q = 


dx 


du 

dy 


Here, we will show how to extend the methods described in Chapter 2 to equations that are nonlinear in p and/or q; that 
is, we consider equations for which f p + fq is a function of p and/or q. (The equation linear in p and q produces an 
expression for f p + f which, of course, depends on only (x,y,u) .) 


3.1 Geometry again 

We start, as before, with the specification of the solution surface in the form 


z = u(x,y ) or u(x,y)-z = 0; 



What if 
you could 
build your 
future and 
create the 
future? 


One generation’s transformation is the next’s status quo. 
In the near future, people may soon think it's strange that 
devices ever had to be “plugged in." To obtain that status, there 

needs to be “The Shift". 
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in this latter form, the normal to this surface is ( U x , U y , — 1) i.e. (p, q, — 1). As we have seen, the quasi-linear equation 
(expressed as ap + bq — C = 0 ) has a tangent plane which is completely determined by the vector (a,b, c) ; this is no 
longer the case for the general equation. Now, at any point on a characteristic line, p and q are related by the requirement 
to satisfy the original PDE: 


f(x,V,U,p,q) = 0. 

To be specific, consider the point (x 0 , Vq , Uq ) in the solution surface and the associated values p — Pq and q = C/q ; 
these five quantities are necessarily related by 


/(*b> v 0 ,w 0 , j Po^o) = 0 ' 

which provides a functional relation between Pq and CJq . 

Sufficiently close to the point X = Xq, y = yn, Z = Uq , the tangent plane in the solution surface takes the form 

z-u 0 =(x-x 0 )p 0 +(y-y 0 )q 0 , 

where we assume that Pq and CJq are not both identically zero. (If they are both identically zero, we have a trivial, 
special case that we do not need to pursue.) Now as Pq and q 0 vary, but always satisfying f (xq , }’q , Uq , Pq , ) = 0 , 

the tangent plane at (Xq ,Vq,Uq), being a one-parameter family, maps out a cone with its vertex at (Xq , JPq , Uq ): this is called 
the Monge cone. [G. Monge, 1746-1818, French mathematician and scientist (predominantly chemistry), made contributions 
to the calculus of variations, the theory of partial differential equations and combinatorics, but is remembered, mainly, for his 
work on infinitesimal, descriptive and analytical geometry. He was, of his time, one of the most wide-ranging scientists, working 
on e.g. chemistry (nitrous acid, iron, steel, water with Lavoisier), diffraction, electrical discharge in gases, capillarity, optics.] 

Some characteristic lines, and associated Monge cones, are depicted in the figure below. 
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Then a specific, possible tangent (solution) plane (which will touch a Monge cone along a generator), in the 
neighbourhood of a particular Monge cone, is shown below. 



3.2 The method of solution 

We are now in a position to describe the construction of the system of ODEs that corresponds to the solution of the 
general, first order PDE 


f(x,y,u,p,q) = 0; 

this method embodies the characteristic structure that we met in Chapter 2, together with the requirements to admit the 
Monge cones. Because we now start with a point on the characteristic lines, but add the requirement to sit on a plane 
through the generator of the cone, the solution will be defined on a surface in the neighbourhood of the characteristic 
line. This surface, which contains a characteristic line, is called a characteristic strip. As we shall see, this process involves a 
subtle and convoluted development that requires careful interpretation, coupled with precise application of the differential 
calculus. The most natural way to proceed is to describe the solution of f = 0 , for p and q, given x, y and u in terms of 
a parameter (as we used in the discussion of the Monge cone): we set p = p(j) , q = q(j) , so that 

f(x,y,u,p(r),q(r)) = 0. 


The T -derivative of this equation yields 


^L^p_ + ^L^L = o 

dp dr dq dr 


which is more conveniently expressed, upon the suppression of the parameter, as 


fp + fq 


dq 

dp 


= 0 


on a particular cone where the PDE provides q = q(p) (since x, y and u are fixed). 


(A) 
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In addition, following our earlier discussion, consider a solution which is described by the solution surface represented by 
a parameter s: X = x(s) , V = y(s) , Z = z(s ) . A tangent vector in this surface is then given by (x^s), y'(s), z'(s )); 
this surface, containing a solution of the PDE, is z = U (x, y) i.e. U (x, y) — Z = 0 , which has the normal (ll x ,u y , — 1 ) 
or, equivalently, (p, q, — 1 ) . Thus we have 


(x',y',z')-(p,q,~ 1) = 0 


which gives 

px' + qy' — z' = 0 or z = px + qy '; 


(B) 


this is an alternative representation of the Monge cone, at a point on a characteristic line, mapped out for some q = q( p) 
. Again, on a particular cone, which uses q = q(p) for fixed (x', y', z') , we take the p derivative of equation (B): 


x' + y'^- = 0. 
dp 


Thus, on a Monge cone, we have 


d q _ x _ dx 
dp y' dy' 



In the past four years we have drilled 

A k 

81,000 km 


That's more than twice around the world. 


Who are we? C 

~. -> 

We are the world's leading oilfield services company. Working ^ 

n J 

globally—often in remote and challenging locations—we invent, 
design, engineer, manufacture, apply, and maintain technology 
to help customers find and produce oil and gas safely. 


Who are we looking for? 

We offer countless opportunities in the following domains: 

■ Engineering, Research, and Operations 

■ Geoscience and Petrotechnical 

■ Commercial and Business 

A > 

If you are a self-motivated graduate looking for a dynamic career, 
apply to join our team. 

What will you be? 

careers.slb.com 

Schlumbergen 
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and when this is combined with equation (A), we obtain 

= _i (o 

d p d y f q 

We return to (B) and now interpret this equation in a different way: first 

, , , d z dy 

z = px +qv becomes — = p + q — 
dx dx 

which holds on a Monge cone with tangent planes on which z = U (so on a solution/tangent plane we have Sz = dll ) i.e. 

(D) 

lx 

Thus, with (C), we have 


dw dv 

— = p + q —— on a cone, 
dx dx 


dw fq r dw .. 

- = r+<, T o,f p - = P f p+q f t . 


(E) 


The calculation thus far has been explicitly associated with a Monge cone, on which solutions z = w(x, y) must sit. But 
solutions clearly must also satisfy the original PDE: 


f(x,y,U, p,q) = 0 ; 


we take 8/dx of this equation: 


du 


dp dq 


fx + fu fp o "E fq ~ 

dx dx dx 


= 0 , 


dq d z u dp 

and then with -=-=- (where we assume that W V1 , exists for our solutions), we obtain 

dx dxdy dy ** 


fx + Pfu = -(Pxfp + Pyfq) ■ 


(F) 


Also, again introducing s (as for (B)), we have 


p 0 ) = 


^-(x(s),y(s)) 

dx 


and so 


ds = “ xx ds +Uxy ds =PxX ' + P} ' y '' (G) 

Now we form p' j ( f x + pf u ) , by combining (F) and (G): 

P' _ PS + Py? 
fx E Pfu Pxfp E Pyfq 
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and then use (C) to eliminate f ' : 


fx + Pfu 


p x x '+p } y 

r dv r 

Pjp+Py^U 


/» 


Px x '+p y y' 


Pxfp + p 


y 

y x' J 


fv 


Px x ' + p y y' 

Pxx' + Py y') 


fv 


A corresponding result, obtained by taking d/dy , gives 


(H) 


fy + qfu = ~(P y fp + <l y fq ) 


= ~((lxfp +( l y fq)' 


cf. equation (F). The development that produces (H) can be followed through to give 

q' y' 

f y +qfu fq 

When we elect to write (E) in the form 


Pfp+qfq fp 

and then combine with (C), (H) and (J), we obtain the set 

r r ! r t 

x _ y _ u _ p _ q 
fp fq Pfp+qfq fx+Pfu fy+qfu' 
or, upon the suppression of the dependence on the parameter s: 


(J) 


dx _ dv _ dw _ dp _ d q 

fp fq Pfp+qfq fx+Pfu fy+qfu ( } 

This is the final, required set of (four coupled) ODEs that represent the solution of the PDE. (These ODEs are sometimes 
referred to as Charpit’s equations, although, strictly, this title implies only when they appear in a different context within 
the solution framework for the PDE; see §3.4.2. They are, more properly, the characteristic equations for the general PDE, 
first derived by Lagrange.) 
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Comment: The case of the quasi-linear equation corresponds to 

/ = a(x,y,u)u x + b(x,y,u)u v - c(x,y,u ) = ap + bq -c, 

and then the general system above becomes 

dx _ dv _ d u _ dp 

a b ap + bq a x p + b x q -c x + p(a u p + b u q -c u ) 


d q 


a yP + by q ~ C y+ b( a U P + V/ “ C u ) 
The first three terms, with ap + bq = C , become 

dx dv du 


the system derived in §2.1, but written in the form given in §2.2. In the next term - the fourth - we use 

8 

— (ap + bq-c) = (a x + a u u x )p + ap x + (b x + b u u x )q + bq x - (c x + c u u x ) = 0, 
ox 


to reduce that term to 


dp 


dp 


a P x + bq x a P x + bp y 
dq 


similarly, the last term can be written as 
Now we have, for example, 


aq x +bq y 

dx dv dp 


a b ap x +bp. 


but on a solution curve we have 


dp d .. .. dv b 

— = — p(x, v(x)) = p x + p— = p x + p --, 
dx dx - dx a 


thus 

dx dp 


a ap + bp } 


, d p 

becomes ap x + bp v = a — = a 

y dx 


r 


b 

Px+Py- 

v a) 


which is an identity. A corresponding calculation applies to the term dqj(aq x +bq v ); thus the last two 

dx dy du 

expressions in the system of ODEs are redundant: we are left with — = =-, as expected. 

a b c 
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We are now in a position to be able formulate and, in principle, solve the system to find a solution of the PDE; for our 
later development, it is convenient to revert to the use of the parameter, s, and so rewrite the system (K) as 


dx dy d.it .■ .■ d p . . dq , . 

— = / p . — = / 9 > — = Pf p +qf q > — = "(/x + /?/«)> ^ = -(/y+«/„)- 


(L) 


Example 7 

Solve the system of ODEs that describe the solution of U x ll y = U i.e. pq = U . 


Here we have f (x, y,u,p, q ) = pq — u = 0 , so that 

fx = fy = 0 , f u = - 1 , fp=q,fq=P ; 


the system (L) of ODEs is therefore 


dx dv dt/ , dp dc/ 

— = q, -f- = p> ~r = 2 p<i ’ — = P’ — = q- 

dy dy dy dy dy 
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The fourth and fifth equations give, respectively, 

p = A C' V and q = Be S (A, B arbitrary constants); 
d U 2s 2s 

the third equation then becomes -= 2ABq~ , and so U = ABq~' + C (where C is another arbitrary constant). The 

ds 

first two equations can now be integrated directly: 

x = Bq s + D and y = Ae s + E , 
where D and E are the last two arbitrary constants. 


However, we are not yet in a position to produce a solution - general or otherwise - without some further considerations. 
We do note, in this case, that we may write 

u = Ae s Be s +C = (y- E)(x -D) + C 

it is then a simple exercise to confirm that this function, u(x,y ) , is indeed a solution of the original PDE, for arbitrary 
D and E, but with C = 0 . 

Comment: The solution described by equations (L) has been developed on the basis, in the main, of a geometrical 
argument (which may appear rather obscure in places). However, now that we have these equations, we can readily confirm 
that, together, they are equivalent to the original PDE. To demonstrate this, consider a solution expressed in terms of a 
parameter s, and then construct 

Q = ^-/OO), v(s),u(s), p(s), q(s )) 

= fx X> + f y y +fu U ' + fpP + fq<l' ■ 

Write u\s ) = — u(x(s), v(s)) = u x x' + u,,v' = px' + qv' 

ds ' 

(which corresponds to the first and second equations of the set, used in the third), then we obtain 

Q = ifx+pfu K+ if y + qfu )y' + f P p'+f q f 

= ifx + pfu ) x ' + ify + qfu )/ - ifx + pfu )fp - ify + Pfu )fq 
when we use the fourth and fifth equations of the set. Finally, we write 

Q = O' - fp )ifx + Pfu ) + (/ - fq )ify + Pfu )f 
= 0 
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by virtue of the first and second equations. The result of using all five equations is therefore 

Q = — = 0 and so / = constant, 
ds 

and the PDE selects this constant to be zero: the set (L) satisfies the PDE. 

It is now clear that we can use the construction described above to derive the set of underlying ODEs directly. This requires 
pairing off terms, and introducing factors, so that dQj cl.V = 0 when appropriate choices are made. The only remaining 
issue in such an approach is to confirm that the choice of ODEs constitutes a consistent set for all the unknowns. 

3.3 The general PDE with Cauchy data 

The task now is to seek a solution of 


f(x,y,u,p,q) = 0 


subject to 


u = u(t) on x = x(t), y = y(t). 


However, if this initial data is to be consistent with the PDE, we require suitable choices to be made for p and q, i.e. 


f(x(t), y(t), u(t), p(t),q(t )) = 0 ; 


(M) 


further, with u(t) = u(x(t), y(t)) , we must have 


dw dx 


— = «„-t -u 

d t dt 


y 


— = px' + qy'. 
dt 


(N) 


Equations (M) and (N) are used to determine p(t) and q(t ), consistent with the given initial data and the PDE. This 
information, in conjunction with the initial data itself, provides all the necessary information for the integration of the 
characteristic equations from, say, 5 = 0. 


Example 8 

2 

Find the solution of U x li y = U , subject to U = t “ on X = t, y = 1 + 1. 


The initial data requires 


? d U _ 

pq =t and — = 2 1 = p + q , 
d t 
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and so pit) = q(t) = t. From the solution developed in Example 7, we see that 


p = te s and q = te s , 


both satisfying the initial data on 5 = 0. Then we have U = t e _s with X = /C' V and V = /C' V +1; thus, eliminating t 

2 2 

and s - which cannot be done uniquely - we obtain U = X" or U = (y — 1 ) or W = x(y — 1 ) , but only the third option 
satisfies the original PDE; hence u(x, y) = x(y — 1 ) . 


2 2 

A second example, based on a classical equation, is provided by U x +U V = 1. This equation - the eikonal equation - 
arises in geometrical optics (and we have normalised the speed of light here). 

Example 9 

2 2 

Find the solution of U x + U~ = 1 , subject to It = Xt on X = y = t , where A, is a constant. 


2 2 

The equation is f = p~ + q~ — 1 = 0, and so we have 


fx=f y =f«=0'fp=2 P ,f q =2q. 
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The initial data must satisfy 


22 , dw . dx dv 

p + q = 1 and — = A = p - \-q — = p + q ; 

df d t d t 


in view of the first equation of this pair, it is convenient to introduce p = sin OC, q = COS OC , and then we must have 
sin OC + COS OC = / 1. Thus solutions exist only if A, is such as to allow the determination of OC (a real parameter). The 
solution is described by the set of ODEs: 

— = 2p, — = 2q, — = 2p 2 +2q 2 = 2, — = — - 0. 
ds ds ds ds ds 

Directly and immediately we see that 

p = sin a, q = cos a, u = 2s + At, x = 2s sin a +1 , y = 2s cos a+t, 
all satisfying the initial data on S = 0 . The last pair of equations yields 


S = - 


x-y 


2(sin«-cos«) 


and t = 


x cos a-y sin a 
cos a - sin a 


(and we shall assume that OC ^ 7t/4,57 t/4 , so avoiding the zero in the denominator); the solution is therefore 

„ „ v-x + A(xcosa- vsina) 

cos a - sin a 
_ x(A cos a - 1) + y{\ - A, sin a) 
cos a - sin a 

_ x(sin a cos a - sin 2 a ) + t(cos 2 a - sin a cos o) 
cos a - sin a 

= xsin« + ycosoc. 


ft is clear that we have developed an accessible method for solving the Cauchy problem for general, first order PDEs. 
Of course, any satisfactory and complete outcome of such an integration process requires that we are able to solve the 
various ODEs that arise - but that is a purely technical detail. There is, however, a more fundamental issue that we must 
now address: in what sense, if any, does this approach produce the general solution? Even if this is the case - and we 
must hope and expect that it does - we do not have a precise definition of the general solution. Furthermore, are there 
any other solutions not accessible from the general solution? (This possibility, e.g. singular solutions, is a familiar one 
in the theory of nonlinear ODEs.) Thus we now consider the solution of the general PDE from a different perspective. 

3.4 The complete integral and the singular solution 

Our starting point is to present a formal definition of the various solutions that, in general, are possessed by the general, 
first order PDE: 
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f(x,y,u,p,q) = 0 . 

Consider the function 

tfr(x,y,u,a,b) = 0, 

which is treated as defining a function u(x, y,Ct,b ) (implicitly), dependent on two independent variables, (x, v) , and 
on two parameters, ( a,b ). We take, separately, the x and y partial derivatives, to give 

^frx frli^X ~ 0 an d (fry + (fru U y — 0 ■ 

These two equations, together with (fr = 0 itself, can be used to eliminate a and b, to produce a first order PDE 

f(x,y,u,p,q) = 0 (p = u x ,q = u y ); 

if this is our PDE, then <fr(x,y,U,a,b ) = 0 is a solution. (The elimination of a and b is possible provided that the 
appropriate Jacobian, <fr xa <fry b — (fry a (fr x p > is not zero, which we assume in this case.) Thus we have a two-parameter family 
of solutions, on the basis of which we can be precise about the nature of all solutions of the PDE. 

1. A solution that involves both independent parameters, ( a,b ), is called the complete integral. (The 
construction of a complete integral will be described in detail below.) 

2. A solution that is described by b = b{a ), for arbitrary functions b, is called a general solution. This can be a 
particular case of the general solution, by choosing a specific b(a) , which can be further particularised by 
choosing a value of a, or by constructing the envelope generated by varying a (see below). 

3. If an envelope of the solutions in (1) exist, then this is a singular solution. (This solution is obtained by 
eliminating a and b between 

(fr a = 0, (fr b = 0 and (fr{x , y, U,a,b) = 0; 

In case (2), the procedure is, of course, adopted for the one parameter a, but this is no longer a singular solution.) 

The fundamental question now is: how do we construct the complete integral? Once we have developed these ideas, we 
shall be in a position to discuss all these possible solutions. 

3.4.1 Compatible equations 

The procedure that leads to the construction of a complete integral requires, first, the notion of compatible equations. We 
are given the PDE that we wish to solve: 


f(x,V,U,p,q) = 0. 

Let there be another equation, involving a parameter a, which has a solution that is also a solution of f = 0 ; we shall 
write this equation as 


g(x,y,u,p,q,a) = 0. 
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The equations f = 0 and g = 0 are said to be compatible. 

These two equations can be solved for p(x,y,U ) and q(x,y,ll ) (provided that J = d(f, g)/d(p,q) ^ 0 ); then 
we have a p and q which are consistent with a solution ll = u(x, y) . A convenient way to express this is to consider a 
solution described by a parameter s: 


u(s) = u(x(s),y(s)), 

du dx dv . , 

and so — = u Y - \-u., — =px+qy. 

d9 y ds ^ ^ 

This can be expressed in the equivalent form of an exact differential as 

du = pdx + qdy 

the integral of this equation generates a second arbitrary constant, b, and hence a solution </>(x,y,ll,a,b) = 0. 
Additionally, in order that pdx + qdy — du = 0 be completely integrable - as it must be - there must exist suitable 
conditions on p and q ; we now find these. 


Let the integral of 


pdx + qdv - du = 0 
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be y/(x,y,u) = 0 ; the differential form of this is 

<// x dx + iff y dy + i// u du = 0, 


and these two versions must be identical. This will be the case if 


y/ x =ap, \ff y = aq, y/ u =-a , 

where Oc(x, V,u) (^ 0) is an arbitrary function. Thus we obtain the identities 

Yxy = ( a P)y = («9)x; V xu = ( <*P) U = ~°x- Vyu = («?)« = ~ a y ■ 

which give the set 

a y P ~ a x q = a(q x - p v ); a u p + a x = -ap u ■. a u q + a y = -aq u . 

Eliminating OC x and OC v yields 

p(-<*q u - v a u ) - q(-ap u - p<* u ) = a (q x - p y ) 

i.e. qp u — pq u + p y —q x = 0 , 


which is a necessary condition that ensures integrability. This equation will provide the basis for constructing a compatible 
equation g(x, y,U, p,q,a) = 0 and then, finally, a complete integral. 

3.4.2 Equation for g 

The given equation f (x, V, M, p(x, y, u), q(x, y, «)) = 0 , when we take 8/8x and, separately, 8/dll , gives the pair 
of equations 


fx + fpPx + fqVx = 0 and fu + fpPu + fcflu = 0 ■ 


Now we form f x + pf u : 


fx + Pfu = -(fpPx + fcflx + PfpPu + PfqVu ) 
= ~(Px + PPu )fp - (Px + pq u )fq ■ 


The corresponding calculation for g(x, V, U, p(x, y, ll), q(x, V, ll)) = 0 gives 


g x + pg u = ~(Px+pp u )gp - (q x + pq u )g q ' 


and then eliminating (p x + pp u ), we obtain an expression for (q x + pq u ) : 


q x + pq u =j[(fx+pfu )g P - (g x + pg u )f P 
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d(f, g ) 

where J =-= f v g„ ~ f c ,g v ^ 0 (as previously assumed). 

d(p,q) 

The calculation is repeated by taking djdy and then d/dU , which produces 

Py + <lPu = j[^Sy + qgu )fq ~ (fy + <lfu )gq 
These two identities are now used in the integrability condition 


Py+qPu = c lx + P c lu 


to give the equation 


(gy + (lg U )fq ~ C fy + Pfu )gq = (fx + Pfu )gp ~ ( gx + Pgu )fp • 


This is the equation for g, which we may write in the form 

fpgx +f q gy+ (Pfp + qfq )g u ~ (fx + Pfu )gp ~ (fy + tfu )gq = 0 


this is a linear equation for g(x, V,U,p,q, a) = 0 . The methods of §2.4 provide the solution expressed by the system 
of ODEs: 


dx _dy _ du _ dp _ d q 


J p J q yJ p “ r HI q J x “ r PJ u J y L U u 
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which are precisely the equations for the solution of f (x, V, U, p,q) = 0 , (L)! These - now - familiar equations, derived 
in this way as the equation for g(x, y, U, p,q,a) = 0, are usually called Charpit’s equations. [Not much is known about 
Paul Charpit de Villecourt, a Frenchman; he is often described as a ‘young mathematician; he died in 1784, probably 
‘young’ - his date of birth is not available. He submitted a paper to the French Royal Academy of Sciences in 1784 - 
note! - but it has never been published; in it he described his method for solving PDEs. His work was known to a very 
small group of mathematicians in France, his methods first being presented by Lacroix in 1814. His work is, ultimately, 
a recasting and alternative interpretation of Lagrange’s.] 

At first sight it might appear that we are no further forward; we have merely, by a devious route, produced the same 
solution method as for our original PDE (f (x, y,U, p,q) = 0). However, there is an important difference here: any 
integral of this system is sufficient to define a compatible equation, g(x, y,U,p,q,a) = 0. Once we have this equation, 
we may solve it with f (x, V, U,p,q) = 0 for p and q; finally we integrate 

du = pdx + qdy 

to find a complete integral of our original equation. (Note that, although any integral that generates a g will suffice, it 
must necessarily contain p and/or q.) 

Example 10 

Find a complete integral of U x ll y = U . 


We have f (x,y,u, p,q) = pq-u = 0 , so that 


fx =fy=°'fu = ~ l Jp ~ d' fq ~ P 

dx _ dy _ du _ dp _ dq 
q P 2pq p q 

The last equality has the integral p/q = a , where a is an arbitrary constant; this is therefore a compatible equation. This 
equation, with pq = U , yields the solution 


P = 




(signs ordered). 


Then 


dw = pdx + qdy becomes dw = ±(y[au dx + ^u/a dy) 


which can be written 


du 



±{y[a dx + —!=dy) or 2\fu = ±(\[a x + — \=y)±2b , 
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where b is the second arbitrary constant of integration (and the additional —2 is merely a convenience). Thus we have 
a complete integral 



where a (> 0) and b are the two parameters in the solution. 


Comment: This solution should be compared with that obtained in Example 7. A general solution is described bv arbitrarily 
assigning b = b{a ) . Both the a and b partial derivatives, each set equal to zero, requires that b + 4- (yfa X + ^=T) = 0 
i.e. U = 0 is a singular solution. That this is a singular solution is confirmed by the observation that there is no choice of 
the two parameters, a and b, that recovers the zero solution. 

Let us explore these various properties of solutions by considering another example. 


Example 11 

2 2 

Find a complete integral, a general solution and a singular solution of U x + lt~ =1 + 2 U . 


We have f = p~ + q" — 2u — 1 = 0 , so that 


f x =f v =0,f u =-2,f=2p,f=2q; 


the system of ODEs is 


dr dv 


d u 


dp d q 


2 p 2 q 2(p 2 +q 2 ) ~2p -2 q 

The last pair here give, directly, p / q = a and then, with the original PDE, we obtain 


7 9 l + 2w \l + 2u 

(1 + a~)q~ =1 + 2 u and so q — -y , p = - - (signs ordered). 


1 + 0 “ 

Thus du = pdx + qdy becomes dll = + 


l + 2u 
'l + o 2 


1 + 0 

(odx + dv) , which gives 


dll 


■ = + 


adx + dy 


\j 1 + 2u yjy 


and so v 1 + 2 u = ± 


ax + y 


+ o" 


a/i 


±b’ 


+ o 


where b is the second arbitrary constant (and the additional + is a convenience). We have a complete integral 


1 

u = — 
2 


r 


ax + y 


\ 




+ b 


1 

2 


v v f + o y 

The general solution is represented by setting b = b(a) . The singular solution (requiring the d/da and d/ db derivatives 
each set to zero) has the property 
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ax + y 

yjl + a 2 


+ b 


= 0 , so that the singular solution is U = 



The only significant issue that we have yet to address is how we apply Cauchy data to a solution described in terms of a 
complete integral. This will also enable us to explain more fully how the designation b = b(a) is used. 

3.4.3 The complete integral with Cauchy data 

Suppose that we have a complete integral, (/>{x, y,U,a,b) = 0 , of the general equation f (x, y,U, p,q) = 0 . Further, 
let the Cauchy data be represented by 


u = u(t) on x = x(t), y = y(t ); 

a solution of f = 0 , satisfying this data, will be one of three possibilities (provided that a solution does indeed exist). 
These are 

1. A particular and specific choice of the values for a and b. 

2. The singular solution (which, of course, is completely determined). 

3. A choice of b{a ), followed by the construction of the envelope solution generated as a varies. (A choice of 
b(a ), followed by a specific value chosen for a is equivalent to case (1).) 

Cases (1) and (2) are exceptional and are easily checked in particular exercises; the case that requires careful analysis is 
(3). This is the situation that arises most often. 

Once we have excluded the possibility that specific values of a and b can admit the given initial data, we have excluded 
particular versions of the complete integral. However, another possible solution will be one in which the data curve touches 
each of the surfaces associated with the complete integral or, at least, a subset of them. The requirement that the Cauchy 
data be consistent with the complete integral gives 

</>(x(t),y(t),u(t),a,b) = 0, 

and that it touches each of these implies that any solution for t must be a repeated (double) root. This is most easily 
expressed by imposing the additional condition 

^-</>(x(t),y(t),u(t),a,b) = 0. 

£it 
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The elimination of t between (j) = 0 and d(f)/ d/ = 0 yields a relation between a and b: this prescribes b = b(ci) . 
However, we have already excluded the relevance of specific values of a and b, thus we must generate another solution 
consistent will all the preceding requirements. This is the envelope generated by the variation of a, given the b(a ) that 
is consistent with the initial data that generates double roots for t. (The situation that corresponds to a single root is when 
the data curve intersects a surface that contains a complete integral; this is the case when specific values of a and b allow 
the initial data to be satisfied.) Geometrically, this means that the initial-data curve everywhere touches all - or a subset 
of - the surfaces that correspond to the family of complete integrals: the envelope of this complete integral. 

We shall explore these various aspects in the next two examples. 

Example 12 

2 

Find the solution of XU x U y + yu~ = 1 which satisfies U =t on X = 2t, y = 0 . 


2 

The equation is f = xpq + yq -1 = 0, and so 

f x =pq.f y =q 2 .fu =°'f P =xq.f q =xp + 2yq ; 
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the system of ODEs is then 


dx _ dy _ d u _ dp _ d q 

xq xp + 2yq xpq + q(xp + 2yq) pq q 2 


The last pair of terms integrate to give p/q = a - a compatible equation - and then the original PDE requires 

2 1 ci 

(i ax + y)q — 1 and so q = + , = , p = + , =■ (signs ordered). 

y]ax + y yjax+v 

Now du = pdx + qdy becomes dii = ± <2( j X + ^ i.e. u = ±2 -yjax + y + b; 

yjax + y 


A complete integral can be written as 


(w —by = 4 (ax + y), 


where a and b are the two parameters. 


It is immediately clear that there is no choice of a and b that permits the initial data to be satisfied by the complete integral, 
so we require the appropriate envelope solution. 

The initial data requires ( t — by = 8 at , and the double root for t that t — b = 4 a (being the f-derivative of the preceding 
equation); the elimination of t yields b = —2a . Thus the relevant complete integral, with this h((l) , becomes 

(.u + 2a) 2 = 4(ax + v), 


and then the envelope solution requires 4 (u + 2d) = 4x i.e. U = (x — ll) / 2 . The solution is therefore given by 


X 2 = 4 


^ X-U^ 


V z 7 


x + 4v or u(x,y) = 


x~+4 y 
2x 


Example 13 

2 2 

Find a complete integral of y(u “ — U “ ) + UU v = 0 , and then the solutions that satisfy 

2 

(a) u=3t on x = 2 1, y = t;{ b) u =2t on x = t “ , y = 0 . 

The PDE is f = y(p~ —q 2 ) + uq = 0 , so that 

f x =0,f y =p 2 -q 2 ,f u =q,f p = 2 vp, f q = -2 yq + u, 

and the system of ODES is therefore 


45 


Download free eBooks at bookboon.com 

















An introduction to partial differential equations 


The general equation 


dx _ dv 


d u 


dp 


d q 


2 yp -2 yq + u 2yp z + q(u - 2yq) pq p 2 -q 2 +q 1 

dp dq 22 

The last pair of terms are -=-, which integrate to give p —q = a , 

q p 

which is a compatible equation. Thus the PDE becomes ay + uq = 0 , so 

ay 


q 


and p =+. a + 


2 2 
a y 


.2 ' 


Now du = pdx + qdy becomes 


J 2 2 

a + Q ' civ - — d v or 
u 2 u 


du + ^dy 

11 = ±dx, 


which is rewritten as 


udu + aydy 

r 2 2 

+ ay 


= ±J~a dx. 


a + 


2 2 
a y 


This can be integrated directly, to give 

sji/ 2 +ay 2 = ±\[a x±b or u 2 = (b + xja) 2 -ay 2 ; 
this is a complete integral, with parameters a and b. 

(a) The initial data, applied directly to this complete integral, yields 

(3f) 2 ={b + 2t4af -at 2 . 

which is satisfied with b = 0 and a — 3; we have a solution expressed as 

2 q s 2 2 \ 

u = 3(x -y ). 

(b) In this case we need a suitable envelope solution; first use the initial data in the complete integral: 

4 1 2 = (b + t 2 yfa) 2 ; 


2 1 

r = 


the f-derivative (for a double root) of this is 8/ = 4 t{b + 1~ yfa)yfa . Thus (selecting t* 0) 

' 2 7 1 

and so b = —j= . 

ya 

The complete integral, with this b(a) , becomes a 1 = (l/ yfa + xja) 2 - ay 2 ; the a- derivative (for an envelope) 
requires 


f 2 ^ 

4 

( 2 ^ 


" 2 ^ 

r b 

W a ) 

,then —j= 

ya 

r b 

W a J 

— 



0 = 2 


1 


+ xV a 


Wo 


1 1 1 

+ - 


2 3 / 2 2 yfa 


-y 
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2 2 1 1 
so y~ = x -— or a = 

a 1 

where the positive root is selected, ensuring that a is real. 



The solution, satisfying the given Cauchy data, is conveniently expressed as 

/ \2 


u 2 = 


1 + 




2 2 

x -y j 


l 


2 2 
X -v 


y 


l 


2 2 
x -y 


This concludes the presentation of the standard techniques for solving the general, first order PDE. Some exercises follow, 
with answers provided at the end of the text, which should help the reader to master these ideas. 

Exercises 3 

1. Find a complete integral of UU x U y = U x + U y . 

2 

2. Find a complete integral of U x +U y = XU . 

2 2 

3. Find a complete integral of x(u x + U~ ) = UU X , and then the solutions satisfying(a) U = 5t on X = 2t , 
y = t; (b) u = 2t on x = 0, y = t 2 . 

4. Find a complete integral of 4 UU X — U y = 0 , and then that solution which satisfies U = At on X = 0 , 

y — t ■ 


** * * * * ** * * * * ** * * 

********** 
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Answers 


Exercise 1 

(a) quasi-linear; (b) general; (c) semi-linear; (d) linear; (e) general. 


Exercises 2 

2 2 

l. u = x + F(y(l-ux)); u=x + exp[j(l-ux)]. 


1 


1 


2 . U = 


F(x~ - e 2y ) - y exp(e z/ - x z -1) - y 


-; u = 


2 y 


3. U =F(y — l/X~ ) ; (a) U = ^ ± —, j ——— — 1 (note: can be written explicitly); (b) does not exist. 


2x x V 4x 


Exercises 3 

l. m = \jlx + a +yj2y + b . 


2. U = ■ 


1 


-(x-a) 3 / 2 ±yfa y + b 


3. u 2 =a x 2 + (ay + b) z ;( a) u = y/5^jx z + y z ; (b) (u z -2y) z = 4(x z +y z )- 


2 , ..2 


.2 , 2 > 


4. W = 


x/a (ax + y) + b 


-|2 


; (a) u =^-(sJy 2 +I2x-y) 
ox 


7 (^y 2 + I2x-y) + y 


+2 


+ 


6x 


yjy 2 +12x- v 


and the boundary condition is satisfied as 
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canonical forms 
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List of Equations 


This is a list of the types of equation, and specific examples, whose solutions are discussed. 

a(x, v) — + b(x, v) — = 0 (general theory).p.53 

ox ' av 

2 

vu x +x u y =0 .p.57 

u x + xu y =0 .p.59 

('j ii f'j if 

a(x, v, u) -1- b(x, v, u) — = c(x, v, u) (general theory).p.59 

ox ' ay 

XU x + UUy = U . . p.61 

d 2 u ? d 2 u 

—— - c —— = 0 (general theory).p.63 

3 r dx z 

withn(x,0) = p(x), u t (x, 0) = q(x) .p.65 

and with u(x,0) = p(x) = 0,u t (x,0) = q(x) = y-e . P-66 


—y - c 2 —y = 0 with u - a sin(ry t ) & u x = 0 on x = Act . 

dt~ 3x“ 

a(x,y)u xx +2b(x,y)u xv + c(x,y)u vy = d(x,y,u,u x ,u y ) (general theory) ... 

u xx - k Uyy = 0 (characteristic lines). 

Uy = ku xx (characteristic lines). 

u xx +Uyy = 0 (characteristic lines) . 

y~u xx -4 x Uyy = 0 (characteristic lines & canonical form). 

2 2 

x"u xx + 2xyu xy +y u yv = 0 (characteristic lines & canonical form). 

2 

y ti xx T 2-XVUxy + |x + 4x ^yy — u x ^ (t +x + 4x j Uy 

(both characteristics & canonical form) 

(l - M 2 ) </> xx + <pyy = 0 (classify). 


p.69 

P-71 

p.74 

p.74 

p.75 

P-77 

p.78 


p.80 

p.84 
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<p x <p xx +0yy = 0 (hodograph & classify). p.85 

[a 2 - u 2 - luvtpyy + [a 2 - v 2 j(/> vv = 0 (hodograph & classify). p.87 

u, + uu x + h x = 0 ; h, + uh x + hu x = h t + ( uh) x = 0 (hodograph & classify). p.89 

a(<p x , <t>y )(f>xx + 2 b{<p x , </>y )</> xy + c(<p x , </> y )</>yy = 0 (general hodograph). p.90 

(' a ~ ~ Qx j xx — 2 (fixfiytpxy ([ a ~ ~ Qy)^yy ~~ ^’ QxQxx ftyy ~ ® > 

u t + uu x + 2cc x = 0; c t + uc x +\cu x = 0 ( all hodograph). p.91 

u t +uu x +2cc x =0; 2c t +2uc x +cu x = 0 (Riemann invariant). p.95 

u t +uu x +2cc x =0; 2c i + 2uc x +cu x = 0 (simple waves) . p.96 
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Preface 


This text is intended to provide an introduction to the techniques that lie behind the classification of second-order partial 
differential equations which, in turn, leads directly to the construction of the canonical form. This topic is likely to be 
included, albeit only in outline, in any mathematics degree programme that involves a discussion of applied mathematical 
methods. The material here has been written to provide a general - but broad - introduction to the relevant ideas, and not 
as a text closely linked to a specific course of study. We start with a brief overview of the relevant ideas associated with 
first-order equations (and there is far more detail in Part 1); this is then developed and applied to second-order equations. 
The main intention is to present the material so that it can be used as an adjunct to a number of different courses or 
modules - or simply to help the reader gain a deeper understanding of these important mathematical ideas. The aim is 
to go beyond the methods and techniques that are usually presented, but all the standard ideas are discussed here (and 
can be accessed through the comprehensive index). 

ft is assumed that the reader has a basic knowledge of, and practical experience in, the methods for solving elementary 
ordinary differential equations. This brief text does not attempt to include any detailed applications of these equations; 
this is properly left to a specific module that might be offered in a conventional applied mathematics or engineering or 
physics programme. However, some important examples of these equations will be included, which relate to specific areas 
of (applied) mathematical interest. 

The approach adopted here is to present some general ideas, which might involve a notation, or a definition, or a 
classification, but most particularly detailed applications of the ideas explained through a number of carefully worked 
examples - we present 20. A small number of exercises, with answers, are also offered, although it must be emphasised 
that this notebook is not designed to be a comprehensive text in the conventional sense. 
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1 Introduction 


The study of partial differential equations, and the construction of their solutions, is at the heart of almost every problem in 
applied mathematics (and many problems in physics and engineering). In this text, we will describe the general character 
of an important class of second order equations - a class that includes the vast majority of equations that are of interest 
and practical relevance. This class will be carefully defined and described, as will the very significant classification of the 
equations. This classification has far-reaching consequences for both the method of solution and the nature of the solution. 
In addition, these ideas lead to a standard representation of the equation - the canonical form - which often, in turn, leads 
to the construction of the general solution. 

We shall start by briefly examining a class of simple first-order partial differential equations in order to introduce the 
concept of a characteristic variable. This is then developed further, thereby providing the basis for both the classification 
and the canonical form of second-order equations. These results are interpreted in the context of the three classical 
equations of mathematical physics, but they are also applied to more general - and important - equations. Additionally, 
some ideas and techniques that relate to methods of solution and interpretation (specifically simple waves and Riemann 
invariants) will be presented. 


1.1 Types of equation 


We first describe the types of equation - both first and second order - that we shall analyse in this text. We consider 
functions, u{x,y) , defined by suitable partial differential equations; in the case of first-order equations, these are 
represented by the relation 


f 


u, x, y. 


3u 3u 

ax'ay 


o. 


We should note, at this stage, that we shall limit our discussion to functions of two variables, although some of the ideas 
go over to higher dimensions; detailed methods of solution (e.g. separation of variables) in two and higher dimensions 
are to be found in the text ‘Partial differential equations: method of separation of variables and similarity and travelling- 
wave solutions’ in The Notebook Series. 


We start with linear, homogeneous equations that contain only derivative terms: 

. . 8u du A 

a(x,y) — + b(x,y) — = 0, 
ox oy 


and then extend the ideas to quasi-linear equations: 


a(x,y,u)^- + b(x,y,u)^- = c(x,y, u) 
ox oy 


i.e. linear in the two derivatives, but otherwise nonlinear and inhomogeneous. The second-order equations that we discuss 
are linear in the highest derivatives, in the form 


a 2 u 


a 2 u 


a 2 u 


a(x, y)—- + 2b( x, y) —— + c( x, y)—- = d 


dx 


dxdy 


dy 2 


>4 y,u. 


au 

dx 


au" 

d y, 
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this is the semi-linear equation. (The equation in which a, b and c also depend on u and its two first partial derivatives is 
the quasi-linear equation, which will not be discussed here, although many of the principles that we develop work equally 
well in this case.) Examples of the three types mentioned above are 


du 

x — + (x + y) 
dx 


du 

dy 


= 0 ; 


du du 

X - Yu -= U ; 

dx Gy 


a 2 u 

dx 2 


i x d 2 u 
+ (x + y)-+ x 


2 a 2 u 


dxdy a/ 


du 

xu H-+ y 

dx y 


r du ^ 2 


dy 


respectively. 
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2 First-order equations 

The fundamental idea that we exploit is best introduced via a couple of elementary functions. First, let us suppose that 
we are given any function, f (x) , then it is obvious (and apparently of no significance) that f (x) is constant whenever 
x is constant; this is usually described by stating that f is constant on lines X = constant ’. Expressed like this, we are 
simply re-interpreting the description in terms of the conventional rectangular Cartesian coordinate system: y = f (x) , 
and then y = f (x) = constant when X = constant . To take this further, let us now suppose that we have a function 
of two variables, f (x, v) , but one that depends on a specific combination of x and y e.g. 

f (x, y) = sin(x- y 2 )- 

This function changes as x and y vary (independently), but it has the property that f = constant on lines 
2 

X — V = constant ; these lines are shown below. 



The function takes (in general different) constant values on each line. It is quite apparent that this interpretation of the 
function provides more information than simply to record that it is some function of the two variables. But we may take 
this still further. 

The conventional Cartesian coordinate axes, and the lines parallel to them, are defined by lines X = constant and 
y = constant , which is usually regarded as sufficient and appropriate for describing functions f (x, v) . However, 
functions such as our example above, sin^X— y^j> are better described by lines X — V = constant (perhaps 
together with either X = constant or y = constant ). These should then be the appropriate lines to use (in place of 
X = constant, y = constant ); these ‘coordinate’ lines are those reproduced in the figure above. The essence of our 
approach to solving partial differential equations is to find these special coordinate lines, usually called characteristic lines. 
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2.1 The linear equation 

Here we consider the equation 

. 8u Ju . 

a(x,y) — + b(x,y) — = 0, 
ox oy 

where the coefficients a and b will be assumed continuous throughout the domain (D) where the solution, u(x,y) , is to 
exist. To proceed, we seek a solution that depends on x (say) and E, = E(x,y) ; the solution is defined in 2-space so, in 
general, we must transform into some corresponding 2-space (that is, using two independent variables). The aim, now, is to 
determine the function %{x, v) so that the equation for u becomes sufficiently simple, allowing it to be integrated; indeed, 
we hope that this results in a solution that depends - essentially - on only one variable (namely, E, ). We note, in passing, 
that if we choose E, = y , then we simply recover the original problem (which does confirm that a transformation exists). 

Let us write, for clarity, 

u(x,y) = U[x,^(x,y)] 


then 

du 8U dE dU du dE dU 

-=-1-and -=-; 

dx dx dx dE, dy dv dE, 
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thus the equation for u (now U) becomes 

a(x, y)(u x +^><11 ^r) + b(x, y)^ y u^ = 0, 

where we have used subscripts to denote the partial derivatives. Now we choose £(x, v) such that a^ x + b <£y = 0 
(and we note that is no more than the original partial differential equation!) which leaves the equation for U as simply 

aU x = 0, 

and so provided a ^ 0 throughout D, then 

U x = 0 or U = F(£) i.e. u(x,y) = F[%(x,y )\, 

where F is an arbitrary function; this constitutes the general solution and confirms that we may, indeed, seek a solution 
that depends on one (specially chosen) variable. (It should be clear that arbitrary constants in the solution of ordinary 
differential equations go over to arbitrary functions in the solution of partial differential equations.) The function £(x, y) 
is determined completely when we consider lines <£, = constant for then 

e , dy e A d y % x b(x,y) 

C r H- C v =0 or -=-=-. 

dx } d x £,y a(x,y) 

Note that, since we now know that we may describe the solution as U = constant on certain lines, we may equally write 
directly that U = constant on lines y' = bfa , without the need to introduce ^ at all. Nevertheless, as we shall see, 
the introduction of characteristic lines is fundamental to any generalisation of this technique. So the integration of the 
ordinary differential equation y' = b(x,y)/a(x, y) yields the characteristic lines C(x,y) = constant (this being the 
arbitrary constant of integration), and then U = b'(c) is the required general solution (which is equivalent, of course, to 
u = constant on lines g(x,y ) = constant ). 

Example 1 

2 

Find the general solution of the partial differential equation yu x + X Uy = 0 . 


The characteristic lines are given by the solution of the ordinary differential equation 

-j ry 1 ^ 

y' = '- —( v ^ 0) and so \y = + constant, 

y 2 J 

2 3 

or g(x,y) = 3y — 2x = constant , which describes the characteristic lines. Thus the general solution is 


u(x,y) = F(3y 2 -2x 3 ), 


where F is an arbitrary function. 
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2 

Comment: We can check this solution directly (at least, if F is a differentiable function), for then U x = — 6x F'(C) 
and Uy = 6 yF\%) , so that 

yu x + x 2 Uy = -6yx~F' + 6 yx 2 F' = 0, 
and observe that this does not require the condition y ^ 0 . 

The next issue that we must address is how the arbitrary function, F, is determined in order to produce - we hope - a 
unique solution of a particular problem. 

2.2 The Cauchy problem 

Let us suppose that there is a line, y = B(x) , Xg < X < X \, along which U = V (x) is given (which might be expressed 
in terms of y or x and y); this describes the boundary condition (given on the boundary curve, y = B(x) ) that should 
determine the unique solution of a first-order partial differential equation. This constitutes the Cauchy problem. [A.-L. 
Cauchy, 1789-1857, the leading French mathematician of his day.] However, two cases can arise and they are significantly 
different. First, if the line y = B(x) is one of the characteristic lines, then V (x) must be constant on this line; otherwise 
the solution does not exist: we would have an inconsistency. If it does exist, the solution is undetermined elsewhere; we 
simply know that u = constant on lines £(x,y) = constant. 

On the other hand, if y = B(x) is not a characteristic line then, wherever it intersects a characteristic line, u takes the 
constant value (along that characteristic line) that is determined by V (x) evaluated at the point of intersection. The 
solution is then defined along y = B(x) , Xq < X < X \, and wherever the corresponding characteristic lines develop in 
the ( x , y) -plane; this is depicted in the figure below. 



The solution exists in the region bounded by the two characteristic lines, labelled C; data is 
given on the line y = B(x ), between x = x 0 , Xj. 


58 


Download free eBooks at bookboon.com 





An introduction to partial differential equations 


First-order equations 


Example 2 

Find the solution of the partial differential equation U x + XU y = 0 which satisfies the boundary condition ll = sin V, 
—oo < v < , on X = 0 . Where is this solution defined? 


d y i 2 

The characteristic lines are given by the solution of — = X i.e. y — — X = constant , and so the general solution is 
u(x, y) = F (y — -1 X 2 j- Now on X = 0 we are given u(0,y ) = F(y) = sin y , so the required solution is 


u(x, y) =sin(y-|x 2 )- 


1 2 

The lines y-yj = constant span the whole plane, and the data is on all of they-axis, so the solution exists throughout 
the (x, y) -plane; this is clear from the figure below. 



All the characteristics intersect the y-axis, and so the whole plane is covered. 


2.3 The quasi-linear equation 

The first-order quasi-linear equation is linear in the first partial derivatives only, so its most general form is 

a{x,y, u) ^ + b(x,y, u) ^ = c(x,y, u). 
ox oy 

Although this equation appears to be much more complicated than the previous one, we may still apply the technique 
developed in §2.1; we introduce £ = <^(x,y) and then seek a solution u(x,y) = U\x,%(x,y)\ , so that we obtain 

a (u x +^xU^) + b^ y U^ = c. 


On lines 
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we have 


a£ x +bi^y = 0 i.e. lines given by 


dy 

dx 


b 

a 


aU x = C ; 


thus we now have a pair of simultaneous ordinary differential equations: 

dy b(x,y,U) dU c(x,y,U ) 

— =- and -=-, 

dx a(x,y,U ) dx a(x,y,U ) 

because they-dependence in the second equation is prescribed by the first. 
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Example 3 

Find the general solution of the partial differential equation XU x + UUy = U . 


In this case, we have the pair of equations 


dv u du u 

— = — = — (jt*0) 

dx x dx x 


so we obtain In w = In x + constant or it = Ax 


dv 

and then —— = A or y = Ax + B , 

dx 


where A and B are arbitrary constants. Thus we may express the solution as 

U 


X 


■ = A is constant on lines V — U = B = constant. 


which can be written as the general solution 


u = xF(y - u), 


where F is an arbitrary function. 


Comments: This solution is, in general, implicit for u(x, v) ; only for special initial data will it be possible to write the 

1 1 / 1 \ 

solution explicitly. For example, given that ll = -y V on X = 1 , then — y = F I — y I i.e. F(z ) = z ; thus 

U = x( V - ll) or li = —. 

1 + X 


Further, we can check that we have obtained the general solution by computing U x and Uy (at least, for differentiable 
F), although some care is required: 


u x F xu x F & ciy _ x(i Uy) F 


and so U x 


F 

1 + xF' 


SC Uy 


xF' 

1 + xF' ’ 


it then follows directly that XU x + UUy = U . 


We now use the idea of characteristic lines to explore the corresponding problems for second-order partial differential 
equations, which is the main concern of this Notebook. However, we shall first explore the special - and simple - case 
afforded by the wave equation. When we have seen how the technique for first-order equations is readily employed 
(Chapter 3) for this equation, we will generalise the approach. 
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Exercises 2 

1. Find the general solutions of these partial differential equations: 

a) 3 u x + 2XUy = X ; (b) u x + 2lly = U ; (c) U t + UU X = 0 . 

2. Find the solutions of these partial differential equations that satisfy the given boundary conditions: 

x , _ 2 

a) U x + Uy — U with ll = X on V = —X; 

Uf + III! x = 0 with U = tanhx on y = 0, and then decide if U x is defined for V/ e[0,oo) and 
Vx e (- 00 , 00 ) . 


x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x- 

x-x-x-x-x-x-x-x-x-x- 
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3 The wave equation 

We motivate the more general analysis of an important class of second-order equations (in two independent variables) 
by considering the equation: 

d 2 u 2 d 2 w . 
c- z —5-= 0. 

dt 2 dx 2 

the classical wave equation. This equation arises in many elementary studies of wave propagation; it describes the amplitude, 
u(x,t ) , of a wave as it propagates in one dimension. (Because we prefer a natural choice of notation for distance (x) and 
time ( t ), these have been used here, rather than the more conventional x and y, although we shall revert to these in our 
analysis of the general equation.) The equation contains a positive constant, c, which will play a significant role in the 
interpretation of the resulting solution. 


3.1 


Connection with first-order equations 


To proceed with our wave equation, we make an important observation: 



fd d Yd 

-h C- 

y dt dx Jv 


__ d_' 
dt dx y 


demonstrating that this second-order differential operator can be factorised. Thus the wave equation can be written as 


fd d 
— + c — 
\dt dx 


Y3u 

Jt dt dx , 


= 0 , 


or, equivalently, 


v3t dx 


3u 

~dt 


■ + 


3x y 


= 0 . 


Clearly a solution of the former equation can be obtained from 


du du 

- c — = 0, 

dt dx 


and from 


du du 
— + c — = 0 
dt dx 


when we use the latter. Thus in the first case we have 


dx 


u = constant on lines — = c 


dt 


i.e. u = constant on x — ct = constant 


so u(x,t) = F(x-ct); 


correspondingly, the second gives 
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u = constant on lines x + ct = constant so u(x, t ) = G(x + ct) . 

All this suggests that we seek a solution of the wave equation in the form u(x,t ) = U(<%, //), where g = X — Ct and 
rj= X + Ct ; thus we obtain 


U x =Uz + and U t =-cUz+ cUyj , 


and then 

u tt _ c u xx = c (u^- 2 U^+U 7777) —C (u^ + 2 U £jj +U tjtj) = — 4 c U^. 

Thus our wave equation 

2 

u tt — c~u xx = 0 becomes = 0 , 
and this latter version of the wave equation has the general solution 

U(£,ij) = F(^) + G(rj) or u(x,t ) = F(x - ct) + G(x + ct), 

where F and G are arbitrary functions; this solution is usually referred to as d’Alembert’s solution of the wave equation. [Jean 
le Rond dAlembert (1717-1783) was a foundling, discovered on the steps of a church in Paris. He had a good upbringing, 
financed by a courtier (who was, it is assumed, his father). He trained first as a lawyer, then spent a year on medical studies, 
before eventually settling on mathematics as a career. He was the first to analyse, in 1746, the motion of a vibrating string, 
introducing both his general solution of the wave equation and the decomposition into harmonic components; he also 
gave a convincing derivation of the wave equation for this phenomenon. It is fair to record that he probably initiated the 
study of partial differential equations. He also made important contributions to dynamics - dAlembert’s principle - to 
complex analysis and to celestial mechanics.] 

We now have two sets of characteristic lines, and each member of each set (in this case) is a straight line, one set with 
slope 1 jc and the other with slope — 1 jc (from t = +x/c + constant ); see the figure below 
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3.2 Initial data 

Suppose that we have Cauchy data on t = 0 (the most common situation for the classical wave equation: initial data ) 
of the form 


u(x, 0) = p(x ) , u t (x,0) = q(x) both for -oo < x < oo, 

although we should note that a special case, of some importance, can be included here: p(x) and Cj(x) both zero 
outside some finite domain. (This is usually described as ‘data on compact support’.) The general solution, expressed in 
terms of F and G, then requires 

F(x) + G(x) = p(x) and -cF'(x) + cG'(x) = q(x ), 


and this second condition can be integrated to give 


G (x) - F (x) 


which then leads to x 

F(x)=|p(x)-^fq(x')dx' 
2 2c J 


1 _ 

c 


x 

Jq(x') dx'' 


G (x) = ^- p( x) +^~ fq(x') dx' 
2 2c J 


Thus, in terms of the given data, the solution becomes 


, . 1 
u(x,t) =-■{ 
2 


p( x-ct) + p( x+ct)+- Jq(x')dx' 


x+ct 


x-ct 


which takes a more readily recognisable form if we assume that there exists a Q{x) such that q(x ) = Q' (x) , for then 


1 


1 


1 


u(x,t) =-\ p( x-ct)--Q(x-ct) ^ + -<^ p{x+ct) +-Q(x + ct) 
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This solution, as we have already noted, comprises two components, one constant along lines X — ct = constant and 
the other along lines X + Ct = constant ; these lines, and hence the solution, cover the half plane t > 0, —co < x < co 
. Observe that, with the identification of t for time and x for distance, then the profile (shape) represented by F(x) (at 
t = 0 ) is a wave component that moves with unchanging form along lines X — Ct = constant , which implies that this 
component moves at a speed c (to the right because C > 0 ); correspondingly, G moves at the speed c to the left. These 
two waves do not interact. Further, if the data is on compact support then, after a finite time, the two waves travelling in 
opposite directions will no longer overlap. 

Now consider a point, (xq , /q ) , in the domain where the solution is defined; one of each of the pair of characteristic 
lines cross at this point. These two lines are 

x-ct = Xq- CtQ (= £ 0 ) and X + Ct = Xq + CtQ (= r/o) 
so they meet the datum line, t = 0 , at X = C() and X = TJq, respectively; see the figure below. 



The triangular region bounded by these two characteristic lines and the segment of the datum line is called the domain 
of dependence: it is the region of the plane in which the solution, u, depends on the data given on the line segment from 
X = to X — //(). A particularly illuminating problem that develops this idea is afforded by the following example, 
especially when we take a limit that results in a concentrated disturbance i.e. a disturbance at a point. 


Example 4 

Find the solution of the wave equation with w(x,0) = p(x) = 0 and U t (x,0) = q(x) = e *where a > 0 is a 
constant, and then examine the form taken by this solution as a —> 0 . (Note that 


jf° e x / a dx+rV x / a dx} = J; 

e x / a ' 

0 + i 

_ e -x / a 

}J-oo Jo J 2 


—oo ^ 

■ 


= 1 - 


so that the area is 1 for all a (> 0 ); as a —> 0 the amplitude (height) increases and the width’ decreases: this is therefore 
a model for the delta function.) 
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The solution is, directly, 


U(Xrt) 


1 

2c 


x+ct 

J q( x') dx' 

x-ct 


1 

f x + ct 

e 

x '/ a dx', 

x> ct 

4ac, 

Jx-ct 

1 

f° e* 

//a dx' + 

1 f x+ct 

4ac, 

Jx-ct 

4ac Jo 

1 

fX + Ct Y 

i e 

//a dx', 

X 

A 

1 

o 

r-t* 

4ac , 

Jx-ct 



- ct < x < ct 


J_ 

4c 

J_ 

4c 

J_ 

4c 


^ e -(x-ct)/a_ e -(x+ct)/a^ x > ct 

( 2 -e (x “ ct) / a - e - (x+ct) / a ), -ct < x < ct 

/ e (x+ct)/a _ e (x-ct)/a\ x < _ ct 


SIMPLY CLEVER 
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and as a —> 0 , this gives (for t > 0) 


U( X, t) —> < 


0, X > Ct 

1/4c on x = ct 
1/2c, -ct<x<ct 
1/4c on x = -ct 


0, x < -ct. 


Thus, for any t > 0, the solution in this limit can be represented diagrammatically as 



and at t = 0 we have the delta function centred at X = 0 . The solution is then non-zero in the sector —ct < X < Ct : 



Comment: This sector, where the solution is non-zero, but which emanates from a point on the axis t = 0 , is called the 
domain of influence i.e. every point in this region is influenced (affected) by the disturbance at the point. 

To complete this initial discussion, we consider one more illuminating example that involves a simple application of the 
wave equation. 
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Example 5 

Find the solution of the classical wave equation which satisfies U = a sin(cyf) and U x = 0 on X = Act (where 
0 < A < 1, a and CO are constants). (This corresponds to a disturbance - an oscillation in time - which is generated at 
a point moving at a speed of Ac in the positive x-direction.) 


The general (d’Alembert) solution is 


u(x, t ) = F(x - ct ) + G(x + ct) 


and we require 


a sin(<»f) = F(Act - ct) + G(Act + ct) and 0 = F'(Act - ct) + G'(Act + ct). 


This latter condition can be integrated once (in t) to give 

1 


F[U-1)ct]+—f— G[U + l)ct]=-. 
c(/L-1) c(/l + l) c 


where A is an arbitrary constant. Upon solving for F and G, these yield 


F [(/l-1)ct] = -la(/l-1) sin(ryt) +^(^ 2 -l) A; 
G[(/l + 1)ct] -^a(A + 1) sin(ryt) -l)A, ; 


F ( a) = -ya(A-l) sin 


aoo 


or G(y3) = la(/l + 1) sin 


Thus the required solution is 


U-1)c 

/3co 

U+ijc 


4 


(T-l)A; 
-i(T-l)A. 


u(x,t) = ^aj(1 + /l)sin 


CO 


c(1 + A.) 


(x + ct) 


(1-/1) sin 


CO 


c(l-A) 


(X-Ct) 


Comment: This solution describes the familiar Doppler effect (shift): the wave moving in the direction of motion (positive 
x) has a higher frequency {G)/( \ — A)) than the wave moving in the opposite direction ( CO/( 1 + A) ). This phenomenon 
was first proposed, on purely physical grounds, by an Austrian physicist C.J. Doppler in 1842; his idea was tested in Utrecht 
in 1845 by moving an open railway carriage, at speed, containing a group of trumpeters who played the same note as 
other trumpeters who stood at the side of the track! 
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Exercises 3 

2 

1. Find the solution of the classical wave equation, — C U xx = 0 , which satisfies 

u(x, 0) = 0, U t (x, 0) = -- y (-00 < X < 00 ). 

1 + X 

2. Find the solution of the signalling problem for the classical wave equation (given in Q.l), for which the boundary and 
initial conditions are 


w(x, 0) = 0 and U t (x, 0) = 0 both in X > 0 ; u(0,t) = f (t) for t > 0. 


(You may find it convenient to introduce the Heaviside step function: 


H (x) = 



x < 0 
x > 0 


in your representation of the solution.) 


**************** 

********** 
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4 The general semi-linear partial 
differential equation in two 
independent variables 

The general equation that we consider here is 

a(x,y)u xx + 2b(x,y)u xy +c(x,y)u yy = d(x,y,u,u x ,u y ) 

where a, b, c and d are given; we have written partial derivatives using the subscript notation. (The appearance of the 
‘2’ here is both a small convenience and a standard convention in these equations.) The basic procedure follows that 
which was so successful for the wave equation, namely, to find a suitable transformation of variables. This will necessitate 
the consideration of three cases, which leads to the essential classification of these equations and then to the standard 
( canonical ) versions of the equation (which might then be susceptible to simple solution-methods). 

4.1 Transformation of variables 

Although we eventually require the solution u(x,y ) , we represent this in the form 

u(x,y) = U[g(x,y), rj{x,y)] 


for suitable choices of the new coordinates 

2j{x, v) = constant, r/(x,y) = constant, 


which replace the conventional Cartesian set: X = constant , y = constant . Thus we have, for example, 

U x — jj ' My ~ ‘ryUz + T]yU jj , 


and then U£ and U,j exist provided that the Jacobian J = 2, x T] y — 2, y Vj x ^0 (and note that the choice 2, = X , 
TJ = y - which is no transformation at all, of course - generates J = 1, so some 2, , 1) certainly do exist). [K.G.J. Jacobi, 
1804-1851, German mathematician, who did much to further the theory of elliptic functions.] However, we also require 
second partial derivatives; for example, expressed as differential operators, we have 


dx' 


a_ 

dx 


' a a ^ 

dz + ' h 3n 


^ a a Y c a a 

a^ drj\ a^ drj) 


and we may choose to use either the first version, or the second, or a mixture of the two. In particular, we elect to use the 
former when we differentiate 2, x and I] x , but the latter when we operate on dj5)2, and d/d rj ; the result is 


d z _ e 8 d 

„ 7 — ?xx ^ dXX n _ 7 

dx 1 d?j dx z 


+ 2 ^xVx 


2 & 
+ >lx ~ 


ocpil ' drj 


.2 ’ 
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2 / 2/2 

there are corresponding results for d / 8x8y and d / dy . Our original equation now takes the form 


AUg +2 BUfr + CU 1V] = D(£, Tj,U,U^U n ). 


where the coefficients on the left-hand side are given by 


A = a£*+2b£ x £ y +c^; B = a£ x /7 x + b(£ x /7 y +£ y 77 x ) + c£ y /7 y ; 

C=a Til+2bTj x Ti y + CT] 2 y , 

and D is a combination of d evaluated according to the transformation and the first-derivative terms that arise from the 
transformation used on the left-hand side. 


2 

The first observation that we make concerns the coefficients A, B and C; in particular, we form B~ — AC (which, as we 
shall see shortly, naturally arises - or a version of it - in what we do later). This gives 


B - AC = 


i2 


a ^x J 7x + b(# x *7y + Bx) + c ^y ; 7y j 

~( a fx +2b£ x £ y + c£ 2 )(a7/ 2 +2b/7 x /; y + c/;J) 


a 2 (d^"d^x) + b 2 {^Tly + Zyrix) 2 -^xZyVxVy 

+C 2 {tfylfy - fyV 2 y) + ac(2^ x ^ y ^ x ^ y - - fyljl) 

[2£x*7x(£x*7y + ^ y rix)-^ 2 xVxny ~ 2^x£y 


+ab 


+bc[2^ y 77 y (^ x ^ y + £ y // x ) - 2^ 2 rj x rjy ~ 2/7y£ x £y 
= b 2 (^x^7y - £y>7x) 2 - ac (^x 7 7y - £y 7 7x) 2 = (b 2 - ac)j 2 


where J = £ x l] v — £> y t] x is the Jacobian introduced above. For the transformation from (x, v) to (£, if) to exist, 

we must have J it 0 , and then the sign of B~ — AC is identical to the sign of b~ — ac (which uses the coefficients 

2 

given in the original equation). Thus, no matter what (valid) transformation we choose to use, the sign of B~ — AC is 
2 2 

controlled by that of b — ac , and this suggests that this property of b^ — ac is fundamental to the construction of a 
solution; the intimate connection with the method of solution will be demonstrated in the next section. 


72 


Download free eBooks at bookboon.com 








An introduction to partial differential equations 


The general semi-linear partial... 


4.2 Characteristic lines and the classification 

Let us address the question ofhow to choose the new coordinates, ^ and // ; lines c(x,y) = constant imply that on them 

= _£x 

dx £ y 

(and correspondingly dy/dx = — Tj x jt] y on lines Jj(x,y) = constant ). With these identifications, we may write 


A = £ 


dy 


- 2b 


■'dy' 
dx j 


+ c 


and C = 1j\ 


dy 
dx j 




2b 


/ dy ) 
dx, 


+ c 


and then both A and C are zero if we elect to use as the definition of the characteristic lines 


f dy 

vdXy 


-2b 


V dx y 


+ c = 0,,.^ = i|b± 

dx 




ac 


Thus if b > ac we have two real families of curves (defined by the solution of the ordinary differential equation) and 
we may identify one family as ^ = constant and the other as lj = constant : we have determined a choice of C and 
T] that simplifies the original equation - it now becomes simply 


2 BU^=D. 


2 2 2 2 
Further, it is clear that we have three cases: > ac , b = ac and b < ac , and we should note that b — ac will, 


in general, vary over the ( x , y)-plane, so there should be no expectation that it will remain single-signed. These three 


cases provide the classification. 
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2 

(a) b — CtC > 0 (hyperbolic) 


This is the most straightforward case, as we have just seen. The characteristic lines, £(x,y) = constant and 
r/{x,y) = constant , are defined by the two (real) solutions of the first-order equation 


dy 

dx 


1 

a 



this is referred to as the hyperbolic case, and the partial differential equation is then said to be of hyperbolic type (a 
terminology that will be explained below). 


Example 6 

2 

Find the characteristic lines of the wave equation U xx — k Uy V =0 ( k > 0 , constant). 


2 2 2 

Here we have a = 1 , b = 0 and C = —k , so that b — ac = k >0 (and so the equation is hyperbolic everywhere); 
thus the characteristic lines are given by 




= ±k i.e. y + kx = constant. 


Comment: This recovers, with a change of notation, the results obtained in Chapter 3. 

2 

(b) b —ac = 0 (parabolic) 

We now have only one solution of the ordinary differential equation, because we have repeated roots; we call this the 
parabolic case. To proceed, we choose one characteristic, ^ say, which is defined by the solution of y' = bja ; the other 
is defined in any suitable way, provided that it is independent of the family %(x, v) = constant i.e. it results in ,J 0 
. Typically, the choice Tj = X is made, although other choices may be convenient for particular equations. 

Example 7 

Find the characteristic lines for the heat conduction (diffusion) equation Uy = ku xx (k > 0 , constant). 


2 

First write the equation as kll xx = U v , then we identify a = k , b = C = 0 which gives /C — ac = 0, so parabolic 
everywhere. Thus y' = 0; we may use ^ = y = constant with JJ = X = constant , which is no transformation at all. 
Thus the original equation is already, as one of parabolic type, written in its simplest form. 
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(c) b~ — ac< 0 (elliptic) 

This case presents us - or so it would appear - with a much more difficult situation: the equation defining the characteristic 
lines is no longer real, so we might hazard that no transformation exists in this case, ft is clear that, because we have the 
identity B 2 - AC = (t> 2 - ac)j 2 < 0. then A and C must have the same sign and cannot be zero; thus we choose to 
define the transformation to produce A = C and B = 0 i.e. 

a - C = a [&- ril) + 2b(£ x £ y - /7 x ^y) + ~ '7y) = 0 

and 

6 =a4/7x+ b (£x'7y+£y'7x) + c £y'7y = °- 

Let us define the complex quantity % = ^ + i T ], then we have 

Xx=Zx+ ir lx and =€y+in y > 


and so 


®X x +2bXxXy + ^X y 

= a (?x- J ll) + 2b (£x£y - My ) + c(4 - 7/y) 

+ 2i [ a £x>7x + b (£x>7y + £y nx) + c ^y ; 7y 

= 0 ; 

thus lines %(x,y) = constant are exactly as before: solutions of 

^ = ^ b ±V b2 -acj = ^ b ±iVac-t) 2 | 

However, the solution of this differential equation is necessarily complex-valued (called the elliptic case), so we write this as 

Z(x, y) = g(x, y) + i r/(x, y) = a + ifi = constant 

where CC + i/3 is a complex constant. Thus the choice of the new coordinates is given by g(x,y) = O' = constant and 
//(X , V ) = /3 = constant (both real!) i.e. we follow the procedure used in the hyperbolic case, but here we apply the 
principle to the real and imaginary parts separately. So there is a transformation, even though the characteristic lines, 
defined by the ordinary differential equation, are certainly not real. 

Example 8 

Find the characteristic lines for Laplaces equation: U xx + W )T = 0 . 
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2 

Here we have a=C= 1 and b = 0, so — ac = — 1<0 i.e. elliptic everywhere, and then y' = +i or 
y + VC = constant . Thus we may choose the transformation C = V and // = X (or vice versa); as in the previous 
example, this is no transformation - the Laplace equation is already in its simplest form. 


The simple results obtained in the last two examples lead naturally to the notion of the canonical form. 

4.3 Canonical form 

The general equation, following a general transformation, is 

A U+ 2BU Zjj + CU 1V j — D 

and then the three cases give 

a) hyperbolic ( A = C = 0 ): 2 BU ^ = D ; 

b) parabolic (e.g. A = 0 , then B~ — AC = 0 => B = 0 ): CUjjjj = D ; 

c) elliptic ( A = C , B — 0 ): A^U ££ +U ^ j = D . 
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These constitute the canonical forms (and so we confirm that both ku VY = U x and U xx + U V y = 0 , Examples 7 and 
8, are already in this form). Here, we use the word canonical’ in the sense of ‘standard’ or ‘accepted’. The terminology 
(hyperbolic, parabolic, elliptic) as applied to the classification of partial differential equations, was introduced in 1889 by 
Paul du Bois-Reymond (1831-1889, French mathematician) because he interpreted the underlying differential equation 


v dx 


2b 


dx ) 


+ C = 0 


as being associated with the algebraic form 



- 2b XV + CX 


2 


= terms linear in x and y. 


2 

Then a = C = 0 gives e.g. XV = constant , the rectangular hyperbola; a = b = 0 gives e.g. X = V , a parabola; 

2 2 2 

a cc c , b = 0 gives e.g. X + k y = constant , an ellipse (and a circle if a = C ). 


The construction of the canonical form, via the appropriate characteristic variables, will be explored in three further 
examples (and then we will briefly examine a few specific and relevant applications of these ideas). 


Example 9 

2 2 

Show that the equation y U xx — 4x Uyy = 0 is of hyperbolic type (for X ^ 0, y ^ 0 ), find the characteristic variables 
and hence write the equation in canonical form. 


2 2 2 

We have ' b — ac' = 4x V > 0 forX ^ 0, V ^ 0,.so every where else the equation is ofhyperb olic type. 1 he characteristic 

lines, where the equation is hyperbolic, are given by the solution of the equation y’ = ±-^4X 2 y 2 j v 2 = + 2x/_y 

(and note that V = 0 must be avoided here, anyway) so that y 2 + 2x 2 = constant; we set C, = y 2 — 2x 2 and 
2 2 " 

77 = v" + 2x , to give 


a 

( d d ) 

, d 

( d d ) 

— = 4x 


and — = 2y 

-1- 

dx 

[dn dfj 

dy 

[K J 


Then we obtain 


d 2 „ 

( d d ) 

„ 9 

f d d \ 

( d d \ 

dx 2 

Jv Kj 

+ 16x 


ydrj d^ 


d 2 , 

(a a ] 

- 9 

( a a "j 

(a a 

a y 2 


+ 4y 

1 a// y 
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the original equation becomes, with w(x, v) = U(dj, rf) , 

4y 2 (u,-U f ) + 16x 2 y 2 (u„-2U^+Ug) 

—8x + —16x y (u^ + 2U^+U jjtj) = 0 


Thus 


64x 2 y 2 U ^ = 4(y 2 - 2x 2 )U ^ - 4( y 2 + 2 x 2 )U ^ 


where we now write y = y (£ + rj) and X = ( Tj — d^) , giving 

which is the canonical form of the equation (because the only second-order derivative is U ^). 


Example 10 

2 2 

Show that the equation X U xx + 2xyu YV + V It vv = 0 is of parabolic type, choose appropriate characteristic variables, 
write the equation in canonical form and hence find the general solution. 


2 2 2 2 

Here we have 'b - ac'=(xy ) —XV =0 , so the equation is parabolic (everywhere). One characteristic line is given 
by the solution of y' = xy / X 2 = y/x (so, technically, we must avoid X = 0) i.e. y jX = constant; thus we introduce 
% = x/y (slightly more convenient than y/x) and choose rj = X , to give 


8 _ 1 8 + 8 8 _ x 8 

dx y dd; 8ij 8y y 2 dd; 


Then we obtain 


a 2 

Q 0 

h 

Q 0 

QJ 

h 

QJ 

a 2 _ i a x 

QJ 

h 

QJ 

a 

ax 2 

l y a^ 1 d?i) 

1 ^ 

1 

- 

1 

1 r"D 
\>^ 

dxdy y 2 3^ y 2 

&■ 

ro 

- 

CD 

>x 



0^ = 2x^ x 2 8 2 

8y 2 ~ y 3 d£ + y 4 df ' 

and so the equation becomes, with u(x,y ) = U(<%, rj) , 
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“7 u & + ^ u $i +u nn 
v y ’ 


( i a 

_ 1 ,, X ,, X ,, 

+ 2xy 

v y y y 


+y 


7 ^ + 7 u « 


= 0 


This simplifies to 


U nn = 0 and so U = F(g) + T]G(g) , 


where F and G are arbitrary functions; thus the general solution is 

u(x,y) = F(x/y ) + xG(x/v). 
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Example 11 

Show that the equation 

2 2 1 

y 2 u xx +2xyu xy +(x 2 +4x 4 )u yy = -^-u x + -(y 2 + * 2 +4x 4 )u y 

(for X ^ 0, y 0) is of elliptic type, find suitable characteristic variables and hence write the equation in canonical form. 


We have _ac' = (xy) 2 — y 2 (x 2 +4x 4 j = —4x 4 y 2 <0 f° r x ^ 0> ^ 0, so elliptic and the characteristic 

lines are given by the solution of the equation 


v ,_ *v±v- 4 *y _ x j- 2 i * 2 

V 2 y y 


Thus \y 2 = \x 2 ±\vc 2 + constant or y 2 -x 2 +$rbc 2 = constant; 

2^23 ^ 3 

2 2 3 

We choose £, = y — X and 1] = X , although we could use just lj = X ; the current choice will produce the simplest 
version of the canonical form - indeed, we could even include the factor 4/3 (and we comment on this later). Thus 

d . d , 2 d d _ d 

— = -2x - f3x - and— = 2v -, 

dx dd; dr] dy d% 

and then 


d z y £ d , 

—- = -2 — + 6x-2x 

dx 2 df dr] 


f d yd )d 


-2x—- + 3x— 


di dqjdf 


— +3x' 


f d j d )d 

v d/]Jd/] 


r 


d .,3^ 


with -= 2y —2x——i-3x z —— .. 

dxdy \ d£ dt])dZ dy 2 S# ' d£ z 


d d A „d A 2 & 
— and —— = 2 f 4 y — 


Thus the original equation, with u(x,y ) = U (<%, ?]) , becomes 


y 2 (-2U p + 6x11 „ + 4 x 2 U ^ - 12x 3 U ^ + 9x 4 ' 


£ + ox\j rj + ^x u^-izx u^ + yx 

+2xy(-4xyll ^ + 6x 2 yU + (x 2 + 4x 4 ^2U ^ + 4y 2 U ^ j 


^-(-2xU^ + 3x 2 U^) + 2(y 2 + x 2 +4x 4 )(2yU^) 


which simplifies to give 
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x 4 y 2 (9U^ + 16U^) = 0 or 9U W + 16U^ = 0- 

This equation is essentially the classical Laplace equation (and therefore the required canonical form); it can be written 

3 

in precisely the conventional form if, for example, we replace T] by (3/4j// i.e. the new // = (4/3)x , which is exactly 
the transformation suggested by the solution of the ordinary differential equation. 


4.4 Initial and boundary conditions 

Any differential equation will normally be provided with additional constraints on the solution: the given boundary and/or 
initial data (as appropriate). Indeed, any physical problem or practical application will almost always have such auxiliary 
conditions. However, what forms these should take in order to produce a well-posed problem for partial differential 
equations is not a trivial investigation. (By ‘well-posed’ we mean the conditions necessary for a unique solution to exist.) 
We have already touched on this aspect for first order equations (§2.5) and for the wave equation (Chapter 3); we will 
now discuss these ideas a little further (although it is beyond the scope of this text to produce any formal proofs of the 
various assertions of uniqueness and existence). 

(a) Hyperbolic equations 

The standard type of data - Cauchy data - is to be given both u and dujdn on some curve, F , which intersects the 
characteristic lines i.e. at no point is T parallel to a characteristic line (so a characteristic line and F do not have a common 
tangent at any point). Here, duj dn iy the normal derivative on T (and this situation is exactly what we encountered 
for the wave equation: ll(x, 0) and - (x,0) were prescribed). Further, it is quite usual to seek solutions that move 

dt 

away (along characteristic lines) from the curve T on one side only (again, as we did in Chapter 3, where we had data 
on t = 0 and we required a solution in / > 0). 

A degenerate version of this problem, closely related to the conventional Cauchy problem, is when u alone is prescribed on 
r , but some additional information (usually u again) is given on one of the characteristic lines ( Cq , say). This is possible 
provided that all the data is consistent with the equations and characteristics, because we have the two sets of information 
needed to determine a solution: one set on the other characteristic lines parallel to Cq and intersecting T , and one set 
on the characteristic lines of the other family, provided that they intersect Cq . This is called a Goursat problem. [E.J.-B. 
Goursat, 1858-1936, French mathematician who made some significant contributions to analysis.] 
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(b) Parabolic equations 

It will be helpful, in this brief overview, to consider the canonical form of the parabolic equation, written with x (distance) 
replacing T] and t (time) replacing E, (see §4.3); the simplest such equation is U xx = U t . The characteristic lines, as we 
have seen, appear as a repeated pair defined by dl/dx = 0 (see Example 7); interpreting this in the form dx/d t , we see 
that propagation on the characteristic lines t = constant is at infinite speed, implying that the whole domain is affected 
instantaneously (although often to an exponentially small degree well away from the initial disturbance). Then we may have 
data on t = 0 (initial data) and, if the solution is defined in the domain t > 0, — oo < x < oo , no further information 
is required (although a boundedness condition may need to be invoked e.g. the solution decays as |x| —> oo ). However, 
more often than not, the region is bounded, usually by one or two lines X = constant , although any pair of curves in 
(x, f)-space will suffice to describe the region where the solution is to exist; see the figure below. 



The solution is in D, bounded by the curves r| , T 2 and x i < x < x~, (on 1 = 0). 

Note, however, that no part of the curves, Fj and T 2 , must be parallel to the characteristic lines i.e. no point of these 
curves must have a slope parallel to the x-axis. 

The data given on the curves, Fj and T 2 , will be either u (the Dirichlet problem) or du/ dll (the Neumann problem) 
or a mixture of the two, each on different sections of Tj and T 2 (the mixed problem). (A linear combination of u and 
dll/dn is also allowed; this is usually called the Robin problem.) [P.G.T. Dirichlet, 1805-1859, German mathematician who 
made important contributions to both analysis and number theory; C.G. Neumann, 1832-1925, German mathematician; 
G. Robin was a French mathematician.] 
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(c) Elliptic equations 

This class of problems is the easiest to describe in terms of boundary conditions. First, initial data has no meaning here, 
for the two variables - see §4.3 - appear symmetrically and there are no real characteristics, so there is no exceptional 
variable such as ‘time’. Indeed, elliptic equations in two variables arise exclusively in two spatial dimensions. Then we 
simply need to prescribe u (Dirichlet) or dll/dfl (Neumann), or a mix of these two, or a linear combination of them 
(Robin) on the boundary of a region, D, in order to define a unique solution throughout D. (Note that, by the very nature 
of Neumann data, the solution in this case will be known only up to an arbitrary constant.) 

Exercises 4 

Write these partial differential equations in canonical form. 

a) Uyy - XU xx = 0 for X > 0 ; 

b) 4 yu xx + 2(V- \)u xy ~u yy = 0; 

c) xyu xx + 4x vu xy + 4x yu yy = u ; 

d) u xx + 2u xy + 2(l-2x+2x 2 )u yy = 0- 


x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x- 

x-x-x-x-x-x-x-x-x-x- 
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5 Three examples from fluid 
mechanics 


The study of fluid mechanics has a long and involved history; it has produced, arguably, the greatest range of applied 
mathematical problems of any branch of study - and it certainly encompasses all the types of partial differential equation 
that we are discussing here. Although other fields of study, e.g. gravitation, elasticity, quantum mechanics, to name but 
three, could be cited, fluid mechanics rather readily provides all the examples that we could ever wish for. Thus the Euler 
equation of inviscid fluid mechanics gives rise to both elliptic and hyperbolic problems, and the Navier-Stokes equation 
is essentially parabolic. We shall describe three simple examples, together with some relevant background, that arise 
as suitable models for inviscid fluid flows; so that we can describe a fair amount of detail, these are problems that are 
hyperbolic, or become hyperbolic under certain conditions. 


5.1 The Tricomi equation 

For steady, compressible flow of a gas, we have 

Du 1 D d 

V • ( pa) = 0 and - = - Vn (where -=-l-U-V ( = u ■ V here)) 

Dtp D t dt 

in the absence of body forces; the usual notation has been used, and we shall assume that the gas is modelled by p = kp 1 ' 
(where k and 1 < y < 2 are constants). The two displayed equations are, respectively, the equation of mass conservation 
and Euler’s equation. Let us consider flow in the ( x , y )-plane which is a perturbation from the uniform state: U = (US)) , 
p = P() and p = Pq , where 17, Pq and Pq are constants. We write the solution as the uniform state plus a perturbation 
in the form 


u = (U+U,v) , p = p Q + r, p = Pq + P, 


and treat u, v, P and r (and their derivatives) as small. The three scalar equations, upon the elimination of P (via p = kf? 
), become approximately 

A)( u x+ v y) + u />x =0; Uu x +-^-y- p x = 0; Uv x +^-/? y = 0 

P 0 P 0 

2 

where y Pq j Pq is usually written as CIq , (Jq being the sound speed associated with the uniform state. Let us further 
assume that the flow is irrotational, so that ( U,V ) = W (j) and then p— U /a q j (where we take all perturbations 
to be zero if one is). Thus we obtain 

Po(</>xx + 0yy ) - (A)U 2 /a 0 ) ^xx = 0 or (l - M 2 ) ^ xx + ^ yy = 0. 

where M = U /<2q is the Mach number of the flow. We see immediately that, for M < 1 (subsonic), the equation for 
(f) is elliptic; on the other hand, for M > 1 (supersonic), the equation is hyperbolic. This is our first observation, and 
it is very significant. [E. Mach, 1838-1916, Austrian philosopher and theoretical physicist; in 1887 he published some 
photographs of projectiles, showing the shock waves.] 


84 


Download free eBooks at bookboon.com 




An introduction to partial differential equations 


Three examples from fluid mechanics 


On the other hand, if M is close to unity (usually called transonic flow), the simple-minded approximation that we have 
used above is not valid; in this special case the appropriate approximation results in the nonlinear equation 

fixfixx & yy ~ 0 ■ 


Now this equation is not of the form 


a{x,y) u xx +b(x,y)Uyy =0, 

so our development that was expounded in Chapter 4 is not applicable. But this equation can be transformed into an 
alternative version that is of the form that we have previously discussed; this is possible by virtue of the hodograph 
transformation. This important idea and procedure is described in §5.4, which provides an appendix to this chapter. In 
the case under discussion here, we form the equation for 

X = X(u,V ) and V = Y(u,V ) (where U = (j) x , V = (fry) 

to replace that for (fr{x, v) ; this produces the equation 

^uu ^ ^vv — 0, 

as described in §5.4. This is the Tricomi equation (F. Tricomi, an Italian mathematician, who published his analywsis of this 

2 

equation in 1922) for which 'b — ac' = — U , so the equation is hyperbolic for U < 0 and elliptic for U> 0. These two 
cases correspond to supersonic and subsonic flow, respectively: the fundamental nature of the problem changes across the 
sonic line (U = 0 ) - on one side there are real characteristics and on the other they are imaginary (i.e. the elliptic case). 
Tricomi introduced this problem as a means for analysing this change of type; however, in the context of compressible flow, 
the situation is rather more complicated e.g. just subsonic flow past an aerofoil produces a region - whose determination 
is part of the problem - where the flow is supersonic, as depicted in the figure below. 


just subsonic flow 

-► 



region of 

supersonic flow 


aerofoil 


(A Tricomi-type problem was encountered in Exercises 4.) 

5.2 General compressible flow 

We suppose that a flow of a gas (described by p = kfX , as in §5.1) satisfies Bernoulli’s equation, in the form 

a 2 +^-(^-l)(u 2 + v 2 ) = constant, 
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where a{p) = ^Jdp/dp = ijyp/p = ^kvp i ' ^ , the local sound speed, and U = (u, v). Further, we assume that the 
flow is irrotational (so U = (j) x , V — (f) v , again as in §5.1) and that the usual mass conservation applies: 

H x +Hy = 0 - 

all written in two-dimensional Cartesian coordinates. Thus we may obtain 

2aa x +(7-l)(uu x +vv x ) = 0 , 


and correspondingly for the y-derivative, where e.g. 


a x = Ur-^^p ir 1 )/ 2 p x 


which gives p x = —-^-(uu x + W X ). 


Hence we may write 


M X + Hy =-^~( UU X +W X )U + /7U X —^-(uu y +W y )v + /7V y 


P 



What if 
you could 
build your 
future and 
create the 
future? 


One generation’s transformation is the next’s status quo. 
In the near future, people may soon think it's strange that 
devices ever had to be “plugged in." To obtain that status, there 

needs to be “The Shift". 


86 


§0 


Click on the ad to read more 


Download free eBooks at bookboon.com 

















An introduction to partial differential equations 


Three examples from fluid mechanics 


and so we have the equation 

(a 2 -u 2 )u x - uv(u y +v x ) + (a 2 -v 2 )v y = 0 

Note that this derivation has not required us to make any approximations about, for example, a uniform state (as we did 
in §5.1): within this model, this is exact (and so, not surprisingly, it is highly nonlinear). 

Now we write U = (f> x , V = (f> y in the derivative terms (only) and so obtain 

(a 2 - u 2 )^ xx - 2 uv 0 xy + (a ' 2 - v 2 )^ yy = 0 ; 

it is convenient to define ® = X(j) x + V(f> v —</> = XU + yv — (f) (usually called a Legendre transformation) and then to 
apply a hodograph transformation: ® = ®(«, V’) (see §5.4). This eventually produces 

(a 2 -u 2 )® w + 2uv® uv + (a 2 -v 2 )® uu = 0, 

for which 

'b 2 -ac'= (uv) 2 - (a 2 -u 2 )(a 2 -v 2 ) = a 2 (u 2 +v 2 - a 2 ), 

so the equation is hyperbolic if i< 2 + V 2 — a 2 > 0 i.e. (u 2 + v 2 )/a 2 = M 2 > 1 (supersonic) and elliptic if 

2 2 2 n 2 V // 

U + V —a <0 i.e. M <1 (subsonic). [A.-M. Legendre, 1752-1833, French mathematician who worked on 

number theory, celestial mechanics and elliptic functions.] 

This example, because it avoids any ‘small disturbance’ approximation, explains completely the hyperbolic/elliptic nature of 
the inviscid, compressible flow problem as it appears in classical fluid mechanics. One quite notable result of this analysis, 
in the hyperbolic case, is that weak disturbances to the flow field, it turns out, propagate along the (real) characteristics 
(usually called Mach lines in this context); larger disturbances produce shock waves. These weak disturbances give rise 
to minute changes in density which, nevertheless, can be measured and photographed: the characteristic lines of the 
hyperbolic problem can be observed! This phenomenon is shown in the photograph reproduced below. 
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This shows a flat wedge in a supersonic flow, Mach number 1 • 8; it has an included angle of 10-2° and its centre line is 
inclined at an angle of 3 • 3° . The faint lines, mainly emanating from the upper and lower walls of the wind tunnel are 
the characteristic lines for this problem/flow; they are visible because of the small blemishes on the surface of the walls. 
Across these lines is a very small change in the density of the air which, in turn, changes the refractive index of the air; 
this change can be captured photographically, by a technique known as Schlieren photography The darker lines are the 
paths of the shock waves. 


5.3 The shallow-water equations 

The assumption of shallow water or, equivalently, long waves, in the classical water-wave problem implies that the pressure 
( p ) in the flow relative to the hydrostatic pressure distribution is proportional to the amplitude of the surface wave. 
Thus, written in non-dimensional variables, we set p = !), where Tj(x,t) is the wave height measured relative to the 
undisturbed surface of the water; thus we obtain 


D u 
D t 


—p x which becomes Uf + UU X + WU Z 


~ r /x 


(all non-dimensional) with the equation of mass conservation for an incompressible fluid: 

u x + w z =0. 


The boundary conditions are 

D/7 

W = = rj t + ut! x on Z = 1 + 1) and W = 0 on z = 0 , 

where the free surface is represented by Z = 1 + p and the impermeable bed is z = 0 . The solution that is relevant here 
is recovered by writing 
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w= -zu x (x, t) = z 


f m + U rj x A 


i +n 

and so we are left with the pair of equations (which are independent of z): 


u t + uu x + h x = 0; hf + uh x + hu x = h t + ( uh) x = 0, 


the shallow-water equations, where h = 1 + rj(x,t) (the local depth). These equations - which are coupled and nonlinear 
- can be analysed in a number of different ways (and the methods described in the next chapter are particularly relevant 
here), but for the present, let us content ourselves with the following observations. 

Although it is possible to eliminate either u or h between this pair of equations (most easily done by first setting U = (f> x 
), this is not very instructive (and is mildly tiresome); indeed, it turns out to be more useful to retain the pair at this 
stage, but to introduce C = , and certainly h > 0 in any realistic problem. (It transpires that this quantity, c, is 

directly associated with the speed of propagation of the surface waves.) With the use of c, our pair of equations becomes 

u t + uu x + 2 cc x = 0 ; c t + uc x + J cu x = 0, 


and then we apply a hodograph transformation (§5.4) in the form 

t = T(u,c), x = X(u,c). 



In the past four years we have drilled 

A k 
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This gives, for example, 


Uf X c /J and C, XJJ, 


where J = X u T c — X c T u ( ^ 0) is the Jacobian; there are corresponding results for the x-derivatives. All this leads to 


X c - uT c + 2 cT u =0; X u - uT u +\cT c = 0 , 

which are linear equations in X and T; to proceed, we eliminate X - the most convenient (by forming X uc ), to produce 



C 


which has the (real) characteristics U ± 2c - we have a hyperbolic problem. This was to be expected, since we had hoped 
to have formulated a problem that self-evidently described wave propagation. 


5.4 Appendix: The hodograph transformation 

The essential idea that underpins the hodograph transformation is surprisingly simple: it is the interchange of dependent 
and independent variables. This is almost a standard procedure for ordinary differential equations; consider, for example, 
the equation 

[xa(y) + b(y)]^= 1, 


which is, in general, nonlinear, non-separable and not of homogeneous type. However, if we treat X = x(y) , then the 
equation becomes 


—— a (y)x =b(y), 

dv 


which is linear in x, so it can be solved completely using standard (and elementary) methods. 


We consider the second-order partial differential equation, of the form 

<*(tx>ty)txx +2b Wx>ty)txy +c (tx>ty)tyy = 0 


first we introduce the Legendre transformation 


< t> = x$x+y<l>y-<l>' 

and treat ® where a = (f) x and f5 — (fry . Thus we have ® = X(X + v/?— (fr , and so 

^a = x + cxx a + (f/ a — {y. a (fr x + y a (fr y | 


= x + ox a + f}y a -(m a + /}y a ) = x 
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and, correspondingly, ® p = V . Then we find 

lx( 0y5 ) = Ix (y ^ = ° S ° = 0 

where CC X = (j) xx and J3 X = (f) v) ,, which produces the expression 

</>xx®aP + < l>xy®pp =0 - 
Similarly —— (® a ) = —— (x) = 0, and then we have 

dy dy 

fixy^aa + 0yy®aP = 0 • 


thus, for ® a p it 0 , we obtain the identities 


® RR ® 

<f>xx = - ( t>xy and &yy = ~~^ £L ^ ) xy ■ 

Q ap ®ap 


We now transform 


® RR ® 

a( t > .xx ~ ~ a ~Z &xy + 2b(f> X y ~ c ~ $xy 

®ap " Q ap 


and hence, for (f) x) , it 0 , our equation 


®$xx T '^-b(j) X y -\-c(j)yy 0 


becomes 


a(a,P)<&pp-2b(a,P)<& a p + c(a,P)<& aa = 0 ; 

this is now in standard form, and so amenable to the methods described earlier. We see that the equation uses OC = <fi x 
and P = </> v as the independent variables; this is the essence of the hodograph transformation (although, as we have 
seen, the general partial differential equation also requires a Legendre transformation). 

In the early applications of this technique, OC and p were velocity components e.g. OC = U = (f> x and pi ~ V ~ (f) ,,. To 
represent a curve in terms of the velocity components, rather than the conventional position vector, is called a hodograph. 
(This terminology comes from the Greek word 060 c; - ‘hodos’ - which means ‘way or ‘road’ and is used to indicate that 
the curve is mapped out other than in the usual fashion via the position vector.) 

1. As a first example, let us consider the general compressible-flow equation (§5.2): 

(a — (ft X _ 20 x (py(pyy T ^3 ~ (j) y^(f) yy ~ 0 , 

which is of the general type described above. Hence we obtain directly 
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(a 2 - u 2 )o vv + 2 uvO uv + (a 2 - v 2 )o uu = 0 ; 
where 0 = XU + yv- <p and (u, v) = ((/) x , )■ 

2. The Tricomi equation (§5.1) 

0X0XX + <! ) yy = ® 

is initially treated in the same manner; the general result, applied in this case, gives immediately 

;/0 +0 =0 

w VV UU J ' 

But with 0 V = V, we can take advantage of the particular form of this equation, and elect to introduce 

y = O v (u,v) = Y(u,v ); 

thus, when we write 

^uuv ~ ® we obtain Y uu + uY vv — 0, 

which simply uses the transformation of the coordinates alone (and note that X = X(u,v) is not needed here). 

3. Finally, we consider the shallow-water equations, §5.3, in the form 

u t + uu x + 2 cc x = 0 ; c t + uc x +\cu x = 0, 

and because we have not developed a single equation (for u or c), we have no natural choice for a Legendre transformation. 
Thus we exploit the interpretation employed for the Tricomi equation by introducing, directly, only a change of variable: 

x = X(u,c), t = T(u,c), 

which is equivalent to working with a single equation in one variable - either wore- as we did for the Tricomi equation. 
Thus we obtain 

^r(X) = |-(x) = 1 so X u u x + X c c x = 1 
ox ox 

8 8 

and —( X ) = — (x) = 0 so X u U t + X c C t = 0; 

8t 8t 

correspondingly, we have 

3 3 

—(T) = —( t ) = 1 so T..u t + T r c t = 1 
dt dt 
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and T ) = ^-(0 = 0 so T u U x + T c C x = 0 . 

ox ox 

The first and fourth of these equations give 

T T 

U x =— and C x = -— , 

X J X J 

where J = X U T C — X C T U ^ 0 is the Jacobian. Similarly, the second and third equations yield 




x tt 


U, = -and Ct — 

1 J 1 J 


Thus our two shallow-water equations can be written as 


X c - uT c + 2 cT x = 0; X u - uT u +\cT c = 0, 


and then we may form 


X uc = f- (u T c -2 cT u )=|( u T u 4 cT c ) 


which gives 4 T uu - T cc = - T c 

c 


because the cross-derivative term ( T uc ) cancels identically. 


Exercise 5 

Equations that model a slender, axisymmetric jet in an inviscid fluid are 

Uf +UU X = R ~R X and +(^uR 2 

where R(x,t ) is its radius and u(x, t ) is the speed ofthe flow in the axial (x) direction. Use the hodograph transformation 
X = X(u,R ), t = T(u,R) , to show that 

Xr -uTr -R 2 T u = 0 and ^RTr -uT u +X u =0. 

Find the equation satisfied by T(u,R) , and classify it. 



x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x- 

********** 


93 


Download free eBooks at bookboon.com 


An introduction to partial differential equations 


Riemann invariants and simple waves 


6 Riemann invariants and simple 
waves 


We conclude by presenting a brief commentary on two important aspects of the solution of coupled, hyperbolic equations, 
even though the properties that we describe are not always available to us. Let us restrict ourselves to two unknowns in 
two independent variables e.g. our equations for shallow-water waves (§5.3): 

u t + uu x + 2 cc x = 0 ; c t + uc x + J cu x = 0. 

If there exists some function of u{x,t ) and c(x,t ) which is constant along characteristic lines, then this is called a 
Riemann invariant. (The invariant is, of course, the function of u and c that is constant.) Indeed, in a problem such as the 
one we have quoted, we might expect a pair of functions that are constant, one associated with each characteristic. We 
will explore the details, for the shallow-water equations, shortly. [G.F.B. Riemann, 1826-1866, German mathematician who 
probably had no equal in laying the foundations for 20 th century mathematics and theoretical physics.] 


A very significant special case, that makes direct use of the Riemann invariants, is the following. A flow field has at least 
two regions, in one of which the solution can be described by a constant state; the Riemann invariant associated with (i.e. 
on characteristic lines emanating from) this region are therefore all the same constant. In this situation, the problem posed 
on the other set of characteristic lines can be solved completely (and quite straightforwardly). The waves that represent 
the solution of such problems are called simple waves; we will give an example below. 
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6.1 Shallow-water equations: Riemann invariants 

It is convenient, first, to write our two equations in the form Ut + Ull x + 2CC X = 0 ; 2 Ct + 2 UC X + CU X = 0 


which are added to give 


and subtracted to give 


(u + 2 c) t + u(u + 2 c) x + c(u + 2c) x = 0, 


(u - 2c) t + u(u - 2c) x - c(u - 2c) x = 0. 


These two equations are therefore 


|- + (u + c) |-)(u + 2c) =0 and j|- + (u-c) ^-K u - 2c ) =° 

at ax at ax 


where we have the characteristic lines defined by 

dx dx 

— = U + C and — = u — C , 

d t At 

respectively; also note the appearance of U + 2c ; cf. §5.3. Thus we have that 

+ dx 

ll + 2c is constant on lines C : — = U + C 


At 


dx 

and U~2c is constant on lines C : — = U — C ; 

At 


the Riemann invariants for this problem are therefore ( U + 2 c) and (ll — 2 c) . It is useful to express the solution in 
the alternative form 

dx 

U + 2c = F(a) , where CC is constant on — = ll + C (i.e. on C ) 

At 

dx _ 

and U — 2c= G(/3) , where (3 is constant on — = U — C (i.e. on C ), 

At 


where F and G are arbitrary functions. We may note that the characteristic lines here describe propagation that is at a 
speed c upstream/downstream (— C / +C ) relative to the speed of the flow (n). 
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6.2 Shallow-water equations: simple waves 


Let us consider the problem of a wave that is moving to the right into a region (X > 0 ) of stationary water (so u = 0) of 
constant depth /zq (so C = Cq = h() here). Then all the C characteristics emanate from this region of undisturbed flow 


(see figure below) 



and so on the characteristic lines C we have 


u-2c = G(/3 ) = -2cq . 

Thus ll — 2c = ~2 cq is constant everywhere (because this is the value taken on all C characteristics, no matter where 
they sit in the ( x , y)-plane); also U + 2c is constant on each C + characteristic. These two conditions imply that both 
u and c are constant on C + lines, so 

dx 

— = U + C can be integrated to give X — \U + C)t = OC 

d t 

and then U + 2c = F(a) = F\x ~(u + c)t ] . 

Let us prescribe the surface (wave) profile at t = 0 - the initial data - by writing 

h(x, 0) = H(x), 


which is to be consistent with h = Iiq in X > 0 . Thus at / = 0 we have 


F(x) = (u + 2c)\ = (4c - 2c 0 )| = 4 jH(x) - 2 
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and then for t > 0, this gives 


i + 2c = F[x - (m + c)t ] 


4 sJh[x - (zv + c)t ] - 2. 


Thus 4c-2c 0 = 4Vh-2^/h^ = 4JH x-(u +Vh)t -2^/F^ 


which requires that 


h( x, t) = H 


x-(u + ^/h(x,t))t 


with 


u( x, t) = 4^ H x - (u + ^)t] - 2^/hb - 2^/h( x, t) 

= A4h-2^hQ-24h 

= 2 ylh(x,t ) -V^o]- 

Finally, when we use this expression for u, we obtain the complete (implicit) solution for h{x,t ): 


h(x, t) = H 


x-(3Vh(x,t) -2^/hbJt 


Thus we have produced a description of the solution, in this example of a simple wave, but this exact solution - of the 
original nonlinear, coupled equations - necessarily is implicit, in general. Indeed, if the initial wave profile is one of 
elevation, then the solution (for h, for example) will ‘break’ in a finite time i.e. the characteristic lines will cross, requiring 
some appropriate shock (discontinuity) condition to be invoked. 
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Answers 


Exercises 2 

l. (a) u{x, v) = F(3v-x 2 ) + |-x 2 ; (b) u(x,y) = e x F(y-2x) ; (c) u(x,t) = F(x-ut) 


2. (a) u(x, y) = 


2(s-jQ 
y 2 -x 2 +4 


„. , , , . , seclr (x - ut ) 

(b) u(x,t) = tanh(x-wf) and n x =-=- 

1 + /seclr (x - ut) 


, which exists for V/ > 0. 


Exercises 3 


1 


1. u(x,t) = —[arctan(x + ct ) - arctan(x - ct )]. 

2c 

2. u( x, t) = f (t — x/c)H (t — x/c). 


Exercises 4 


(a) 2u 


#7 


- u j^>Z = y- 2Vx , rj = y + 2y[x , hyperbolic; 
2 

(b) 2(1 + = —Ujj, 4 = 2y — x,rj = y + x, hyperbolic; 

(c) rj^+ rj^Ufirj = U + 2r/^4+ f l 2 )' J 4> 4 : => ; -X 2 , I] = X, parabolic; 

(d) + U jjjj - 


2 | 2 

U r j,g = y-x,?i = x -x,elliptic. 


1 + 4/7 


Exercise 5 


27^ m + R^Tjrm + 3R Tft = 0 so 'h 2 —ac'— —2R? < 0 i.e. elliptic. 
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List of Equations 


This is a list of the types of equation, and specific examples, whose solutions are discussed. 


Method of separation of variables: 

^tt — C ^xx ~ ® . 

iif — ku xx . 

u xx u yy ~ 0 . 

— ^(O^xx . 

u t = k[ii n . +ju r \ . 

u rr T7 u r + Uffff — 0 . 

' r 

U rr + 7 U r + -T- + -i- (cot = 0 
w r — k(uxx +^yy) . 

^// ~ C ) _ 0 . 

Urr ^ f Uf "f 2 " Uffff + — 0 . 

/ r 


u f — k 


( 1 ^ 

u /r + “ + “V M 0<z> + “T ( cot + T- 7 u 69 

r r sin (j) j 


r r 1 2 

' r 


n= 


n YX — u tt = sm u . 

Travelling-wave solutions: 


Uft u xx 0 


u i — ku xx 


u t = ku xx in x <ct, with 11 = 0 and u x = —c both on x = ct 


Uxx Uyy ~ 0 


11 XX “f Uyy + U zz 


Uff C ^ 

U t — k(u XX “t Uyy ) 


) =0 


p.107 

p.109 

p.109 

p.110 

p.110 

p.112 

p.112 

p.114 

p.115 

p.116 

p.117 

p.118 

p.119 


p.123 

p.123 

p.124 

p.124 

p.125 

p.126 
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u t + uu x = ku xx . 

u t + 6 uu x + u xxx = 0 



Similarity solutions: 

Uj — kll xx . 

Utf ~ C ^xx ~~ ^ . 

^ xx + ^yy — ® . 

u t + u xxx = 0 . 

Uf — k(u xx +Uyy ^ . 

u t + uu x = ku xx . 

u t + 6 uu x + u xxx =0 . 

u t = (uu x ) x which satisfies J u(x,t) dx = constant 

u t = ku xx (more general similarity solution). 

u t + u xxx = 0 which satisfies u(x, 0) = f(x) . 


p.127 

p.128 

p.130 


p.132 

p.134 

p.136 

p.139 

p.139 

p.141 

p.142 

p.143 

p.144 

p.146 
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Preface 


This text is intended to provide an introduction to the standard, elementary methods that are used for the solution of 
partial differential equations (mainly of the second order). Some of these methods and ideas are likely to be mentioned 
in a degree programme that includes mathematical methods, probably first met in the second year. The material has been 
written to provide a general - but broad - introduction to the relevant ideas, and not as a text closely linked to a specific 
course of study. Indeed, the intention is to present the material so that it can be used as an adjunct to a number of different 
courses or modules - or simply to help the reader gain a deeper understanding of these important techniques. The aim is 
to go beyond the methods that are presented in a conventional study, but all the standard ideas are discussed here (and 
can be accessed through the comprehensive index). 

It is assumed that the reader has a basic knowledge of, and practical experience in, the methods for solving elementary 
ordinary differential equations. In addition, it is assumed that the reader is also familiar with eigenvalue problems; these 
are not developed here, although the ideas will be mentioned. This brief text does not attempt to include any detailed 
applications of these equations; this is properly left to a specific module that might be offered in a conventional applied 
mathematics or engineering or physics programme. However, some important examples of these equations will be included, 
which relate to specific areas of (applied) mathematical interest. The use of transform methods is not included. 

The approach that we adopt is to present some general ideas, which might involve a notation, or a definition, or a method 
of solution, but most particularly detailed applications of the ideas explained through a number of carefully worked 
examples - we present 32. A small number of exercises, with answers, are also offered, although it must be emphasised 
that this notebook is not designed to be a comprehensive text in the conventional sense. 
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1 Introduction 


The formulation of virtually every problem in modern - and classical - applied mathematics results in partial differential 
equations (of various types). The ultimate aim is then, usually, to obtain explicit solutions relevant to the investigation being 
undertaken. Of course, quite often, no solutions can be found (even assuming that appropriate solutions exist) because 
the equation, or the set of equations, are too difficult to solve. Although this may be the eventual outcome, the starting 
point for any such considerations is the collection of standard techniques that enable simple equations to be solved. In 
this volume in The Notebook Series, we shall give an overview of these methods. It is not our intention to produce, for 
example, all the details of Sturm-Liouville theory that are required for the method of separation of variables; this can be 
found in the volume in this series entitled An introduction to Sturm-Liouville theory’. Rather, we shall outline how the 
basic method unfolds when applied to various fairly common partial differential equations of applied mathematics and 
mathematical physics, and particularly when expressed in various coordinate systems. 


The emphasis here will be on how partial differential equations are reduced to one or more ordinary differential equations, 
by a suitable transformation or assumption about the form of the solution. Of course, quite often, other - more general 
- solutions may exist, obtained by applying less restrictive assumptions; some of these ideas can be found in the volume 
‘Partial differential equations: classification and canonical forms’. In this volume, we will develop the techniques that 
produce ordinary differential equations by invoking the method of separation of variables, or by seeking solutions of 
travelling-wave or similarity form. As we shall find, one great advantage over more general methods (e.g. as applied to 
general, semi-linear equations in two independent variables) is that solutions can often be found for more complicated 
variants of such equations. Thus we may be able to cope with highly nonlinear equations, equations of higher order or 
equations in more than two independent variables. Indeed, we shall present examples taken from all these categories, 
many of them being important equations in their own right. 

1.1 The Laplacian and coordinate systems 

The simplest class of second-order partial differential equations that we commonly meet in classical applied mathematics 

2 

and theoretical physics is based on the Laplacian ( V ■ V = V ) of the unknown function. So we have 

V 2 U = 0, U t = kV 2 U and Uff -c 2 V 2 u = 0, 


which are Laplace’s equation, the heat conduction (diffusion) equation and the wave equation, respectively; here t is time 
and k and c are positive constants, and we have denoted partial derivatives by subscripts. The physical systems that might 
be modelled by any one of these equations could be written in any suitable coordinate system; typically, these will be 
either rectangular, Cartesian coordinates (x, y, z), or cylindrical coordinates ( r, 6, Z ), or spherical coordinates ( r, 6, (j) 
), although other choices are possible, perhaps even general curvi-linear coordinates. The simplest, and most familiar, is 
rectangular, Cartesian coordinates, giving rise to the operator 


V 


2 


a 2 a 2 a 2 

ax 2 + dy 2 + az 2 
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Then, with this choice of coordinates and a suitable reduction in the number of spatial dimensions, examples of our three 
equations above are 


8 2 u 8 2 u 

8x 2 + 8y 2 


8u 

8t 




8 2 u 

8x 2 


= 0 , 


which are the simplest equations that are classified as elliptic, parabolic and hyperbolic, respectively. (The background 
and relevant details of this classification can be found in any introductory text on partial differential equations, or in the 
volume in this series: ‘Partial differential equations: classification and canonical forms’.) 
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The other two, standard, coordinate systems are defined by 



for cylindrical coordinates (where X = r COS# , V = /' Sin 6*), and 



for spherical coordinates (x = r sin (j) COS 0 , y = r sin (f) sin 0 , z = r COS (j) ). Note that we have elected to use the 
symbols r and 0 in each, although r has a different meaning in each; there should be no confusion, because either one 
or the other will be used exclusively in a particular application. The corresponding Laplace operators are then 

2 _ S 2 10 IS 2 d 2 

V = -yd-—+ ^r- — H-— 

dr 2 r dr r 2 d6 2 dz 2 

in cylindrical coordinates (which is equivalently plane polars plus z), and 

w2 8 2 2 d Id 2 1 J 1 d 2 

V =—yd—~ + ^- + + ^-y 

dr 2 r dr r l d(j) 2 r 1 r 2 sin^ (j) d6 2 

in spherical coordinates. (We quote these results here; a derivation can be found in any good text on the vector calculus 
which includes a description of general curvi-linear coordinates.) 
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1.2 Overview of the methods 

The reduction of partial differential equations to ordinary differential equations can be accomplished in one of three ways, 
although sometimes we are simply generating the same solution in a different way; it is these three approaches that we 
shall discuss in this volume. The first, and most familiar, is separation of variables; suppose that U = u{x,t ) is defined 
by a partial differential equation, then we seek a solution u(x,t) = X(x)T(t) - the variables have been separated. Now 
this may appear to be a very special construct - it is! - yet quite general problems, based on suitable linear equations, can 
be solved when we invoke, in addition, the method of Fourier series. As we shall see, the method can be extended to any 
number of variables, so we may write e.g. ll(x,y,Z,t ) = X(x)Y(y)Z(z)T(t) . 

The second and third methods are essentially two different variants of just one: instead of separating the variables, all the 
variables are combined together into one functional group. However, there are two very different ways of doing this. In 
one, we restrict ourselves to a linear combination of the variables e.g. (kx + it); this is called the travelling-wave form, 
because it mirrors precisely the type of solution obtained in wave problems e.g. with solutions such as F(x — Ct) , written 
in the usual notation, which is a travelling wave (of unchanging form). (This terminology is quite often used even if time 
is not involved, although any resulting solutions are usually of no great importance.) In the other combination, all the 
independent variables are multiplied together, allowing for general powers e.g. Xt ; this is called the similarity (or self¬ 
similar) form. (Note that F(x m t ^ ) can be replaced by G(xt n ) , for suitable G and n, so there is no need to raise both 
x and t to general powers, in the case of two independent variables.) The solution expressed as F(xt n ) has the property 
that F = constant on lines Xt' 1 = constant; this description of a solution i.e. constant on a given family of curves, 
is usually called self-similarity. A self-similar object is one which essentially looks the same on different scales; a classic 
example is a fractal. In our case, F(xt n ) looks the same for any /'’(constant) , and this applies for any combination 
of x and t (different scales’) for which xt” = constant. 

We should comment that, although the method of separation of variables enables very general solutions to be constructed 
(satisfying rather general types of initial and boundary conditions), the other two solutions are not so accommodating. 
The travelling-wave solution is a special solution of the equation which would normally be generated only by special 
initial (and possibly boundary) conditions, although such solutions might appear eventually (e.g. as time increases) in 
certain cases. This is typically the situation that is encountered in the appearance of solitary-wave solutions from general 
initial data in a ‘soliton problem. The situation that obtains for similarity solutions is even more extreme. Now the initial 
and boundary conditions must be very special indeed, and quite often cannot be easily predicted before the solution is 
known. However, it does sometimes happen that such solutions, suitably interpreted, can be used to derive more general 
solutions of the same equation, but satisfying more general initial and/or boundary conditions. 

Finally we mention that another important method for finding solutions of partial differential equations, also usually by 
generating ordinary differential equations, is to employ transforms. This is too large a topic - there are many different 
transforms - for incorporation in this text; it is hoped to devote a separate volume to a description of transform methods. 
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2 The method of separation of 
variables 


The history of this method is hidden in the fog that inevitably surrounds the past; it was used by a number of mathematicians, 
but usually applied only to very specific problems, from about the middle of the 18 th century Certainly L’Hospital, 
d’Alembert, Daniel Bernoulli and Euler employed the technique, and a case can be made that L’Hospital was the first. 
However, it was not until J.B.J. Fourier (1768-1830) that a complete and systematic development was presented (in 
about 1807, for the problem of heat conduction); thereafter it became a standard procedure in the armoury of applied 
mathematicians, at least for certain types of partial differential equation. 

2.1 Introducing the method 

We will describe the fundamental principles that underpin the method of separation of variables by considering the 
classical wave equation 


U tt — c^u xx =0 ( C> 0 , constant), 

which is of hyperbolic type. The separable solution is written as u(x,t) = X(x)T(t) , for suitable functions X and T ; 
the wave equation then becomes 


XT"-c 2 X"T = 0, 


where the primes denote derivatives with respect to the corresponding arguments of the functions. It is convenient to 
divide throughout by XT: 


T" 2 X" n X" 1 T" 

- C -= 0 or -= —- 

T X X c 2 T 

and, although this manoeuvre requires XT ^ 0, we can dispense with this restriction when we have seen how the method 
proceeds. (We may note that this structure of the equation has produced another version of the separable property.) Thus 
we have 


X" 1 T" 

- (a function of only x) = —-(a function of only t) = A , say, 

X c z T 

but x and t are, by definition, independent variables i.e. they are assigned arbitrarily on the appropriate domains. It is 
clear, therefore, that any one of the choices 

A is a function of only x; A is a function of only t; 2 is a function of x and t, 
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leads to an inconsistency. The only possible choice is A = constant (usually called the separation constant), and so we 
obtain 


X" , T" 2 „ 

X T 

The problem for X(x) now becomes, with suitable boundary conditions, an eigenvalue ( Sturm-Liouville ) problem, which 
will have appropriate solutions only if 

A = -a> 2 < 0 

i.e. X'' + (0 2 X = 0 and then T"+C 2 C0 2 T= 0, 

for specific values of A , the eigenvalues. The general solution for both X(x) and T(t) involve trigonometric functions. 
Finally we observe, from the original separation of variables, namely 

XT" — C^X'T = 0 , that we must have e.g. X" OC X , 

and so T" oc T : the equation is separable without the need to require XT ^ 0 . A more general solution can now be 
obtained, in the familiar way, by summing over all the eigenvalues of the underlying Sturm-Liouville problem (permitted 
because the partial differential equation is linear). 

This method goes over directly to the other two elementary, standard equations of applied mathematics/theoretical physics. 
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Example 1 

Apply the method of separation of variables to the heat conduction (diffusion) equation U t = ku xx (k > 0, constant). 


Weset u{x,t ) = X (x)T(t) .which gives XT' = kX"T and we require for separability and consistency that X" oc X ; 

2 

let X" = XX . The Sturm-Liouville problem, with suitable boundary conditions, then requires that X = —CO < 0 i.e. 

2 2 
X" + CO X = 0 , which (again) has trigonometric solutions. This leaves T' = —kco T , so that 

T(t ) = A exp [-kccrt^, 

where A is an arbitrary constant: the x-dependence is oscillatory, but the f-dependence is a decaying exponential (because 
CO ^ 0 and is real, and k > 0) . 


Example 2 

Apply the method of separation of variables to the Laplace equation U xx + Uy V = 0 . 


We write ll(x,y ) = X(x)Y(y) , and so obtain X"Y + XY" = 0, and then with X" oc X (or, indeed, Y” oc Y), 
we have 


X”+XX = 0 and Y"-XY = 0. 


Depending on the boundary conditions, either X(x) or Y (y) will be described by a Sturm-Liouville problem. If this is 
2 

X(x),then X — CO >0 and X(x) is a trigonometric function, but then Y(y) is a hyperbolic (exponential) function; 
on the other hand, if Y(y) is the Sturm-Liouville problem, the roles of X(x) and Y(y ) are reversed. However, there 
will always be one of this pair of functions that is trigonometric and the other hyperbolic. 


Comment: In summary, we see that we have 

2 

U t f —C It xx = 0 (hyperbolic type); both X(x) and T(t ) are trigonometric; 
ku xx = Uf (parabolic type); X(x) is trigonometric and T(t ) is exponential; 
U xx + Uyy =0 (elliptic type); X(x) is trigonometric/hyperbolic, Y(y ) is 

hyperbolic/trigonometric. 
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The three examples that we have discussed above have been second order, in two independent variables, and with constant 
coefficients. It should be plain that linear, constant coefficient equations are all susceptible to the method of separation 
of variables. However, if the equation is nonlinear, it is not clear if we can proceed (but, exceptionally, this is possible, as 
we shall demonstrate later); if the equation has variable coefficients then, in general, we cannot separate the variables. 
Consider the equation 

u tt -[c(x,t)\ u xx = 0, 


which is a wave equation with variable speed, then separability clearly fails i.e. XT" = c(x,t) X"T is not separable 
unless the c(x,t ) is itself separable, but then we will certainly generate a more difficult pair of ordinary differential 
equations to solve. One example, of this general type, that is readily solved, is the heat conduction equation with a variable 
conductivity: 


U t = k(t)u xx , 


for then we obtain XT' = k{t)X"T , and so X" = -AX yields 

t 

V = -Ak(t)T i.e. lnlZ’l - -Aj k(t') df' 


t 

-A^k(t') dt' . 

0 

rt 

In this case, all that has happened is that the variable t has been replaced by k(t') d t' (provided that this integral 

rt *0 

exist; otherwise we simply work with k{t ') df ' )■ It is instructive to note that we may perform this transformation 
(replacing t by the integral) in the original equation, no matter what method of solution we might adopt. In some of what 
follows, variable coefficients do arise, but by virtue of the forms of the Laplacian that appear in the equations. 


or T{t ) = ,4 exp 


2.2 Two independent variables: other coordinate systems 

We introduce this extension of the method by first considering the heat conduction (diffusion) equation with cylindrical 
symmetry: 

u t = k(u rr + ^ u r j (k > 0, constant). 

The separable solution takes the form u{r, t ) = R(r)T(t) , which gives 

RT' = k(R" + ^R'Y 

and so with T' = —AT , we obtain 

-Aft = £(fr 
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Almost all physically realistic solutions will require that u decays exponentially in time (t), so we shall impose 
A > 0 . We have 


R" + -R' + -R = 0 
r k 


in which we may write R(r ) = S^r^jk/k j , to gi 


or 


give 

-S"+-]£s' + -S = 0 

k r\k k 

pS" + S' + pS = 0 (where p — r^jA/k ), 


which is a zero-order Bessel equation (with solutions Jq and Yq ; see e.g. the volume ‘The series solution of second order, 
ordinary differential equations and special functions’ in The Notebook Series). The solution that is bounded as p —> 0 is 
Jq , so one solution of the partial differential equation is 

u(r,t ) = A J 0 (r-yjA/ k ) exp(-/U), 

where A is an arbitrary constant. (Typically, the eigenvalues are associated with the zeros of Jq or Jq .) 


The same technique can be applied to similar equations, as we now demonstrate. 
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Example 3 


Apply the method of separation of variables to U rr 
polar coordinates). 


+ \, u r + \ 
' r 


Uqq = 0 (which is Laplaces equation written in plane 


We set u(r, 0) = R(r)&(0) , to give 

+ + =0, 

and we write 0” = — 20 ; if the problem is described by 0 < 0 < 2 7t, then we must have solutions that are periodic 

(period at least 2 n ) if they are to be continuous for 9 G [0, 2?r\ . On the other hand, if 0 G [0,@o].0 <0q <2k , 

then no such condition can be imposed; nevertheless, we normally expect A > 0 so that the solutions are trigonometric 

2 

in 9 (rather than exponential or linear). Let us write A = CO > 0, then we obtain 

0(0) = A sin (co0) + B cos (<x>9), 


where A and B are arbitrary constants. 


The equation for R(r ) becomes 


r 2 R" + rR' -co 2 R = 0 


which is of Euler type; the solution, for CO > 0, is of the form R{r) = r a , so 


2 2 2 

a(a-l) + a-co =a -co = 0 i.e. a = ±co , 


which leads to the general solution 


R(r) = Cr a + Dr~ a , 


where C and D are arbitrary constants. Thus we have a solution for u : 


0) = (Cr a + Dr~ eo j(A sin (co0) + B cos (co9 )). 


Comment: We can simplify the dependence on the arbitrary constants; for example, if all the constants in the above 
solution are non-zero, then (for new arbitrary constants) we may write 

u(r,0) = a[i- co + Br~ co j(sin(<y6 r ) + Ccos(co0)). 

Example 4 

2 1 1 

Apply the method of separation of variables to U rr + — U r + —y + — (cot (f^U* = 0 , which is Laplaces equation 

r r y * 

with symmetry about the ‘polar axis in spherical coordinates i.e. u does not depend on 9 . 
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We seek a solution in the form u(r, (/) = R(r)(\R(/j ), to obtain 

(iT + fR'"jd> + ^(RO" + RO' cot </>) = 0, 

2 1 

in which we write R" + — R' = X —r- R , leaving the equation for ® : 
r r z 

®" + ®'cot^ + A® = 0. 

Now this latter equation is more usually expressed in terms of ®(0) = ®(C0S (f)) , so that we have 

®' = -(sin0) v F'( x ) where X = COS0, 

and then 0" = -(cos^'F' + (sin (j)) 2 = -JtH" + (l - X 2 ; 

thus the equation for ® becomes 

(l-x 2 ) v F"-2x v F' + i v F = 0. 

This is Legendre’s equation (see e.g. the volume ‘The series solution of second order, ordinary differential equations 
and special functions’) and, for — 1 < X < 1 i.e. 0 < (f) < In, bounded solutions exist only for X = n(n + 1) (with 
n = 0,1, 2,...); the corresponding solutions are the Legendre polynomials, P n (x) = P n (COS (j )) . 

With this choice of X , the equation for R(r) is 

r 2 R" + 2rR' -n(n + l)R = 0, 

which is again an Euler equation. Thus we seek a solution R(r) = l~ a , which gives 

2 

a(a -1) + 2 a- n(n +1) = a + a - n(n +1) = 0 i.e. a = n or -(n +1) , 
which is valid for all the allowed n (so no repeated roots); the general solution is therefore 

R(r) = Ar n H. 

r n +1 

A solution of the partial differential equation is thus 

u{r,(f>) = ^ Ar n + Br *” +1 ^i^(cos0) . 
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Comment: If r = 0 is in the domain of the solution, then we shall presumably require B = 0 ; correspondingly, if r —> <x> 

, then we shall need to select A = 0 (unless n = 0 ) for a solution bounded at infinity. 

2.3 Linear equations in more than two independent variables 

The method of separation of variables can be extended to any number of variables, at least if the equation is linear and 
with constant coefficients; to see this, let us consider, as an introductory example, the heat conduction (diffusion) equation 
in two spatial variables: 


u t = k[u xx + iiyy j (k > 0, constant). 

First, we write u(x,y,t ) = T(t)(f>(x,y ) , to give 

T'<f) = kT[</) xx +(/> yy ), 

and the separation of the t- and the ( x , y)- dependences requires that T' = —AT (where A is our familiar separation 
constant). This leaves 


k(^xx + $yy ) 
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which is an equation of elliptic type, solved in any appropriate way. But one way is clearly by separation of variables, 
obtained by writing <j){x,y ) = X(x)Y(y) : 

k(X"Y+XY") + XXY = 0, 

and then we may set, for example, X" = jjX (where JJ is a second separation constant) which gives 

kY" + {X + k/Y)Y=0. 

It is now clear that we may seek a solution, a priori, in the form 

u(x,y,t) = X(x)Y(y)T(t); 


further, it follows that we may extend this procedure to any number of independent variables (if the equation admits such 
a solution). Let us return to a consideration of the details of the current equation. 

We expect that X > 0 , so that the solution remains bounded as t —> oo , and then we have two eigenvalue (Sturm- 
Liouville) problems: 

2 2 
first X" + CO X = 0 (where we have written jU = —CO ) 


and then Y + 


(X ^ 

- CO 


F=o,, 

which becomes a second conventional eigenvalue problem if ( Xjk ) — CO > 0 . For example, the boundary conditions 

X(0) = X{a) = 0 and 7(0) = Y(b) = 0 

require the choices CO = nnja (n = 1,2,...) and (X/k) — CO~ = ( nut/by (m = 1, 2, ...) so that we have 

(2 2 \ 


X = kX 


n m 

^" + 72 
y a b j 


>0. 


/ \ 

/ \ 



( 2 

2) 


nn 

— X 

sin 

mn 

xry 

exp 

- kn ' 2 

n 

o 

m 

t 

U J 


Kb ; 



(a 

b 2 ) 



In this case we have a general solution, obtained by summing over all m and n , that can be written as 

00 00 

u(x,y,t) = ^ X A mn sin 

n=lm=l 

where the A mn is a set of arbitrary constants; note that this solution comprises a double Fourier series (in x andy). 
The corresponding wave equation provides us with a closely similar example, as we show below. 


Example 5 


Apply the method of separation of variables to U tl — C~ (u xx + U vv j = 0 (c > 0 , constant). 
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We set u(x,y,t) = X(x)Y(y)T(t ) , which gives 

XYT”-c 2 (X"Y + XY")T = 0, 

2 2 

and then choose to write T" = —CO C T , to leave 

co 2 XY+X"Y+XY” = 0. 

2 

Now we set X" = —Q X , which produces the equation for Y: 

Y" + (ar +Q 2 )r = 0. 

Thus we can have a solution that is oscillatory (i.e. trigonometric) in t, in x and in y. 


Comment: This solution would be relevant, for example, to the vibrations of a rectangular membrane. 

Now we turn to the consideration of three independent variables as they appear in cylindrical coordinates. 

Example 6 

Apply the method of separation of variables to the Laplace equation 

u rr +“ u r + \ u 00 + u zz = ®- 

' r 


We seek a solution of the form u(r, 0,z) = R(r)®(0)Z(z) , then 

R"QZ + ±R'QZ + \Re"Z + RQZ" = 0 ; 

r ,7 

2 

we choose to write 0" = — CO 0 and Z" = AZ , which leaves 

R" + ^R'-\co 2 R + AR = 0. 

r r 2 

This equation, written in the more usual form, becomes 

r 2 R" + rR' + (/ir 2 -ty 2 )R =0,, 

which is a Bessel equation. 
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2 

Comment: Although we are likely to require > 0 (so that Q(0) is a trigonometric function), we might have A < 0 
(so Z(z) is also trigonometric) or A > 0, in which case Z(z) will be an exponential function. These choices then 
control the particular type (and order) of the solution of the Bessel equation (which is described in the volume ‘The series 
solution of second order, ordinary differential equations and special functions’). 


Example 7 


Apply the method of separation of variables to the heat conduction (diffusion) equation 

( . \ 


Ut; = k 


1 

u rr + f U r +^U # + ^(cot^+^-2 

r r r sin 


0 


u ee 


( k > 0, constant). 


We seek a solution in the form ll{r,6,(j),t ) = R(r)Q(0) ( l>((/>)T(t) , which gives 


R0®T' = k 


R"0OT +-R'0®T +4tR0O'T +^^R0OT + „ 1 „ R0"OT 


r 2 sin 2 (p 


further, we set T' = —XT ( X > 0 ) and 0" = —CO~Q , so we obtain 


J 


-ZRO = k 

At this stage, we elect to write 

which leads to the equation for R: 


f 2 \ 

R"0 + -R'0 + ^yR0" + ^^R0'— J° - RO 
r r 2 r 2 r 2 sin 2 0 


co 


®" + (cot^)®'- -z —® = -//®, 

sin" (f) 


-XR = k 


f 2 u ^ 

R" + -R'--^-R 
v r r 2 


The equation for ®(^) = 'F(COS^) is (cf. Example 4) 


(l-x 2 )'f"'-2x'F' + 


CO 


2 \ 


jU 

V 1 - x J 


= 0 ( x = cos (/)), 


which is a generalised Legendre equation; the equation for R(r) is 


r 2 R " + 2r R' + 


a 2 ^ 

-r -ju 

vk 

-V2, 


R =0, 


which can be recast as a Bessel equation (by writing R = r ' S(r ) )■ 
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Comment: Note that the equation for R(r ) is not a Bessel equation, because of the factor 2 r ; to be a Bessel equation 
this must be simply r. The proposed transformation ensures that S(r) does satisfy a Bessel equation. 

2.4 Nonlinear equations 

We conclude our presentation of various applications of the method of separation of variables by showing that, occasionally, 
nonlinear equations can have separable solutions. One particularly simple type of problem, in two independent variables 
( x and t, say), occurs when we have 


L t («) = M x (u), 

where L t is a suitable (even nonlinear) differential operator in t, and M x is a corresponding operator in x. 

Example 8 

Apply the method of separation of variables to ^ j = U xx , p > 0 (constant), which is a nonlinear diffusion equation 

(of a type that can arise in the study of plasmas). 


We seek a solution u(x,t ) = X(x)T(t) , which gives 
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and then we set T' = AT" p i.e. 

fP~ 1 

-= At + A (z?^l), 

P~ 1 

where A is an arbitrary constant. (Note that the case p = 1 recovers the linear equation.) The equation for X (x) becomes 

X" = AX P 

which can be integrated once, by first multiplying by X' : 

X "X' = AX P X ' so X ,2 = ^-( X P+1 + B) 

p+1' ’ 

where 5 is a second arbitrary constant. This allows us to represent the solution for X in an implicit form: 

x r dx' Via 

-J(X ') p+1 + B VP + 1 ' 

which requires A > 0 for a real solution. 


Comment: This example demonstrates that the technique can be used, but we note that the resulting equations are no 
longer likely to be trivial (or even standard). 

A rather more sophisticated type of example arises when we must first transform the equation before seeking a separable 
solution. An interesting (but rather involved) such example is presented in the next exercise. 

Example 9 

Transform the sine-Gordon equation, U xx ~U tl = sin it , according to u(x, t ) = 4arctan[v(x,f)] , and then seek a 
separable solution for v(x,f) (which is most conveniently expressed as v(x,t) = X(x)/T(t )). 


We are given u(x, t ) = 4 arctan[v(x, t )], so we obtain 

4v 


1 + V 


% and U xx 


4v 


xx 


8vv y 


1+v" 


( 1+ v 2 ) 2 ' 
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and correspondingly for the t- derivatives. Thus the sine-Gordon equation can be written 


4(v X x ~ 

1 + V 


v tt) 

2 


8v( 


2 .. 2 ) 

- = sinu = 2 sin(u/2)cos(u/2) 


v x -v t 


H 2 ) 


: 4 sin(u/4) cos( i / 4) cos 2 (u/4) - sin 2 (u/4) 
4cos 4 (u/4)tan(u/4) 1-tan 2 (u/4) 
4tan(u/4) 1 - tan 2 (u/4) 


sec 2 (u/4) 


i2 


4 tan(u/4) 1 - ta n 2 (u/4) 4v(l - v 2 ) 

2 2 

1+tan 2 (u/4) ( 1 + N/2 ) 


Thus we find the equation for v(x,t ) to be 

(v X x-v t t)(l+v 2 )-2v(vJ-v t 2 ) = v(l-v 2 ) = v(l + v 2 )-2v 3 „ 

or (l + V 2 )(v xx - v tt - v) - 2v(v 2 - v t 2 - v 2 ) = 0, 

which does not look promising! However, we do note that, although the original contained both linear (u xx — U^) and 
highly nonlinear (si n u) terms, the equation for v has a more uniform structure. We set v(x, t ) = X(x)/T(t) , to give 


1 + 

f X 1 

2' 

( X" 

1— 
X 

2 XT' 2 

X 1 

2X 


r x'} 

It J 


T 

V 

T 2 

T 3 

T J 

T 


It J 


( XT' ^ 


VT 2 ) 


f Y \ 2 


x_ 

T) 


= o. 


This is more conveniently written as 


2<_ 

T 3 



+ — 

T ) 




= 0 


and then collecting appropriate terms, this becomes 


(X 2 +T 2 )^—2X' 2 + X 2 +(X 2 +T 2 )^-2T' 2 -T 2 =0, 
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which still does not look very hopeful. Although more general solutions are possible, we will content ourselves with a 
simple choice: X' = XX ; let us investigate what this produces. 


First, we observe that X" = AX' = A X ; now, because of the near-symmetry between X and T, we also write T' = juT 
(and then T" = //“ T ) to give 

(l - A 2 +// 2 )x 2 +(a 2 -// 2 -i)t 2 =0. 

2 2 

Thus we do indeed have a separable solution, provided that fl = A — 1; this gives 


X' = AX and T 




or 


X = Ae^ and T=Be 


— IT ( for |/l| > 1), 

±t^h 2 -1 


where A and B are arbitrary constants. A solution of the sine-Gordon equation is therefore 


u(x, t) =4arctan 


C 



where C = AjB and A (with |/l| > 1) are two arbitrary constants (so we do not have a conventional eigenvalue in 
this problem). 


Comment: The solution that we have constructed is a wave, propagating with a speed y/ $ — 1 fA \, either to the left 
or to the right. This is a solitary-wave (so-called ‘kink’) solution of the sine-Gordon equation, one of the important 
‘soliton’equations. 


Exercises 2 

1. Apply the method of separation of variables to the wave equation 

u tt -c 2 (u rr + 7 " u r ) = 0 . 


2. Seek a solution of Utt-c ^u rr 
equation for v. 


+ 2 


r u r ) ~ 0 in the form u{r,t) = r ''’(/%/), and comment on the resulting 


3. Apply the method of separation of variables to U t + UU X = 0 . 


*************** 

X-X-X-X-X-X-X-X-X-X- 
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3 Travelling-wave solutions 

Suppose that we have a function of two variables, u(x,t ) , where it will be convenient to use the familiar interpretation: 
x is distance and t is time. This will enable us to give a description that is consistent with the notion of a travelling wave. 
Further, let us suppose that the function takes the form 

ll(x,t ) = F(x-Vt) , 


for some constant V and some function F, either or both of which may be arbitrary (but it is not unusual to find that they 
are prescribed). We now have a representation that can be expressed as 


dx 

U = constant on lines X-Vt = constant i.e. on lines — = V ; 

d t 


thus all points on the profile (shape), ll = F(x) , move at the same speed (V). This is therefore a wave that propagates at 
speed V (to the right if V > 0), with unchanging shape. (We may note that this is precisely the form of each of the two 
components that constitute d’Alemberfs solution of the wave equation: 


u(x,t ) = F(x — Ct) + G(x + Ct ) , the general solution of U tt 


-c 2 u 


XX 


= 0 .) 
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A solution that can be written as 

u(x,t) = F(x-Vt) 

is called a travelling-wave solution. 

3.1 The classical, elementary partial differential equations 

We have just mentioned the connection with the classical wave equation; indeed, we may use this approach to recover 
dAlembert’s solution. We seek a solution u(x,t ) = F(x — Vt) of 

2 

lift — C U xx = 0 (where O 0 is constant); 

thus we obtain 

V 2 F”-c 2 F" = 0, 

where the prime denotes the derivative with respect to (x — Vt) . Thus for arbitrary F"(x — Vt) ( A 0), we require 

V 2 = C 2 so that V = +C; 

then, because the wave equation is linear, we may use this result to construct the general solution of the wave equation 

u(x, t) = F(x - ct) + G(x + ct) 

where F and G are arbitrary functions (although this argument is only valid for functions whose second derivatives exist 
and are non-zero). 

Example 10 

Seek a travelling-wave solution of the heat conduction (diffusion) equation U t = kll xx (for k > 0, constant). 


We write u{x,t) = F(x — Vt) , where V is an unknown constant, and then we obtain 

-VF' = kF" so that kF' = -VF + A V , 

where A is an arbitrary constant. The equation for F has the general solution 

where 5 is a second arbitrary constant; the travelling-wave solution, in this case, is necessarily an exponential function. 
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Comment: A more usual (and useful) form of this solution is obtained when we set —Vjk = i l (where i is real), 
which gives 

F = A + B exp[if(x + i kit)] = A + B exp^ifx - ki 2 t 

which represents a non-propagating harmonic wave (Cxp( ifx) ) that decays in time. We note that, because the equation 
is linear, we may elect to use this complex form, or take the real or the imaginary part. 

A problem, based on the equation of heat conduction, that more closely corresponds to a travelling wave is described in 
the next example. 

Example 11 

Find a travelling-wave solution of the problem: U t = ku xx in X < Ct ,with U — 0 and U x = — C both on X = Ct , where 
c is a constant. (This is a single phase, one dimensional Stefan problem, which might model a moving front, with e.g. water 
on one side, and ice on the other side of the front X = Ct , and u the temperature (taken to be zero in the ice phase).) 


We seek a solution of the form u(x,t) = F(x — Vt), then we obtain 

—VF' = kF" with the general solution (see above) F = A + Bq ^ /k)( x ^) 

This solution satisfies the two conditions on X = cl if V = C , A + B = 0 and B(V j k) = C , which produces 

u(x,t) = k(e- (c / k)(x - ct) -l). 


Comment: If C > 0 , then in the region behind the front, we have U ^ cc as ( X — Ct ) —> —oc (so the water temperature is 
above zero - usually regarded as the stable configuration). On the other hand, if C < 0 , then U —> —k : the water is colder 
than the ice, which suggests that we have a problem! (It turns out that the existence of appropriate solutions is in doubt.) 

Example 12 

Seek a ‘travelling-wave solution of Laplace’s equation: U xx + Uy V = 0 . 


We write u(x,y) = F(ax + by ) , to give 


a 2 F" +b 2 F" = 0, 
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and so for a and b real (and at least one non-zero), we require F" = 0 i.e. 

u{x,y ) = F = A + B(ax + by ), 

where each of A, B, a and b are arbitrary constants. This is not a very illuminating or important result! 


Comment: If we allowed ourselves the liberty of writing b = ±i a (a real), then we have a solution for arbitrary functions, 
F. Thus we may then write down the complex-valued version of the general solution (cf. d’Alembert’s solution of the wave 
equation mentioned earlier in this paragraph): 

u(x,v) = F(x + i y) + G(x - iy). 

The further choice G = 0 then leads to the fundamental notion of a function of a complex variable, the Cauchy-Riemann 
relations, and so on. 

3.2 Equations in higher dimensions 

The obvious exemplar under this heading is the wave equation written in rectangular, Cartesian coordinates: 

Utt~ c ( u xx +u yy + u zz)-°-- 

The travelling-wave solution then takes the form if(x, V,Z, t) = F\kx + ly + 171Z — Vt ) , 
where k, l, m and V are constants; this gives 

V 2 -c 2 (k 2 + l 2 +m 2 ) 

and so for F" A 0 (but otherwise arbitrary) we require 

V = ±cylk 2 +/ 2 +m 2 . 

This solution describes a wave, of arbitrary shape, that propagates with speed V (either forwards or backwards); the wave 
profile is such that 



u = F = constant on lines kx + Iv + mz-Vt = constant, 

which are parallel planes in ( x , y, z)-space, at any time, t. Such a solution is often referred to as a plane wave. Indeed, if 
we recast the solution as 


u(x,t) = F(k-x-Vt), 


where k = (k, /, in) is the wave-number vector, then lines k • X — Vt = constant can be described by 
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k • — = F = ±cV k 2 +1 2 + m 2 i.e. k • —— = ±c , 
d t d t 

where k is the unit vector in the k -direction. Thus the component of the velocity of the wave, in the direction of the 
wave-number vector, is c (either forwards or backwards), and elementary considerations confirm that the wave propagates 
in the k -direction, with k normal to the wave fronts (defined by the lines k • X — Vt = constant). 

Example 13 

Findthetravelling-wavesolutionoftheheatconduction(diffusion)equationintwospatialdimensions: Uj- = + Uyy ) 

( k > 0, constant). 


We seek a solution u(x, y,t) = F(px + qv — Vt) , where p and q are constants, to give 

-VF' = k(p 2 +q 2 )F", 

which therefore repeats Example 10, but with k replaced by k| p 2 + C] 2 j. 
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3.3 Nonlinear equations 

A class of equations that often pose significant difficulties if we are aiming to find general solutions are nonlinear equations, 
but many of them are susceptible to more direct methods. Here, we will introduce this aspect of travelling-wave solutions 
by considering a classical example of this type: 

u t + uu x = ku xx ( k > 0, constant), 

the Burgers equation. This is an important equation, particularly relevant to models of sound waves in a viscous gas and 
to the structure of shock waves, although it was first introduced as a model for one-dimensional turbulence (by J.M. 
Burgers in 1948). 

We seek a solution u(x,t) = F(x — Vt), and then we obtain 

-VF' + FF' = kF" or -VF + \F 2 =kF' + A, 

where A is an arbitrary constant. The solution usually of interest is the one that satisfies conditions such as 

11 = F —> Uq as X —> co ; u= F —> U\ as X —» -00 , 

which are typically the initial conditions for the unsteady problem i.e. U = Uq in X > 0 and U = U\ in X < 0 , at t = 0 
. These conditions require that the two roots of 

i_ p- — yp — A = 0 are F = Wq , U \, 

so we have F 2 — 2VF — 2 A = (F — Uq)(F — u-|) 

V =1(110 + 11 ,) and A = -lu 0 u q , 


and this expression for V is a significant result: the wave propagates at a speed that is the average of the conditions ahead 
and behind. Now we have the equation 


and so we obtain 


1 

( u 0“Ul) 



where g = X - Vt ; thus 


2kF' = (F-u 0 )(F- Ul ) 


H2 


dF=J 


2k' 


In 


^-“0 

F -u\ 


U °~ Ul 4+B 

2k 


or 


F -u 


F -u-| 


^ = C exp 


f u o ~ °i A 

2k 
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where the arbitrary constant C is most naturally associated with an arbitrary origin shift ( Xq ). Thus we may write 


u 0 + U| exp 


u( x, t) = F(x -VT) 


/ u Q —Ul ^ 

2k 


(x-Vt-Xo) 


1 + exp 


2k 


(x-Vt-Xo) 


which is more conveniently expressed as 


u(x, t) =|(u 0 + u 1 )-l(u 0 -u 1 )tanh 


4k 


(x-Vt-Xo) 


where V — — (uq + U-|an example is reproduced below. (This is often referred to as the Taylor shock profile.) 



It is instructive to observe that, as k —> 0 + (equivalently decreasing the viscosity in the fluids context), so the wave front 
steepens; in the limit we shall obtain a jump (discontinuous) solution, depicted below. 


We record the significant fact that the Burgers equation can be solved completely, by using a transformation (the Hopf- 
Cole transform) that maps the equation into the classical, linear heat conduction equation. It is then possible to find the 
solution, for example, to the general initial-value problem (but this goes well beyond the material of this Notebook). 

Example 14 

Find a travelling-wave solution, satisfying U —> 0 as |x| —> oo , of the Korteweg-de Vries (KdV) equation: 
u t + 6uu x + u xxx = 0. 
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We seek a solution u(x,t ) = F(x — Vt ), then we obtain 

-VF' + 6 FF' + F'" = 0 or -VF -3 F 2 +F" = A; 

we assume that U —> 0 (at infinity) sufficiently smoothly so that U x —> 0 and U xx —> 0 there, and so we must have 
the arbitrary constant A = 0 . Now we form 

-VFF' + 3F 2 F'+F"F' = 0 which gives -^VF 2 + F 3 +]-(F ') 2 = B , 

and by the same argument as above, the second arbitrary constant ( B ) must also be zero. This final equation for F can 

2 

be integrated directly, although this is mildly tiresome; it is simpler to confirm that a suitable SCCll - function is the 
relevant solution. Let us write 


F(£) = a sech 2 ( b£) (g= x-VT-xq) 


where a, b and Xq are constants; we then have 

F' = -lab sech 2 (b<^) tanh(&£), 
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and so our equation gives 


1 

2 


Fa 2 sech 4 (6£) + a 3 sech 6 (6£) + 4a 2 6 2 sech 4 (b<%) 1 - sech 2 (jb%) 


= 0 . 


We therefore have a solution when we make the choices 


V = 4b 2 and a = 2b 2 , 


for arbitrary b (and also arbitrary Xq , which simply represents an origin shift); thus the solution is 


u(x, t) = 2b 2 sech 2 


b|x-4b 2 t- 



Comment: This is the famous solitary-wave solution of the Korteweg-de Vries equation; it plays an important role in water 
waves and most particularly, in ‘soliton’ theory 

We conclude with one more example, of some practical significance. 

Example 15 

Find a travelling-wave solution, U = F(x — Vt), of the equation Ut = ( u n u x ) ( n > 0, (n > 0 , constant), which is zero 
for X > Vt (i.e. ahead of the wave front). 


We write u(x,t) = F(x — Vt) , to give 

-VF' = (F n F') andso-VF = F n F'+A, 

and for F = 0 (and F' finite) on E, = X — Vt = 0 , we require the arbitrary constant A = 0 . Thus we have 


F n ~ 1 F' = -V and so -F n =-V%+B, 
n 

and the same condition just used now requires B = 0 ; we have the solution 

[nV (Vt - x)] 1/n , x <Vt 


1 


u(x,t) 


0, x > Vt. 


Comment: This solution, although continuous, does not have continuous first or higher derivatives if fl > 1, so we would 
need to reconsider the applicability of the equation in this case. The equation describes the (one dimensional) flow of a 
perfect gas through a porous medium, where u is the speed of the gas; the equation, with the given condition ahead, is 
regarded as a good model for 0 < n < 1. 
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Exercises 3 

1. Find a travelling-wave solution of U t = kll xx + Au , where k > 0 and A are constants. What condition ensures that 
oscillatory solutions exist? 

2. Find the equation satisfied by the travelling-wave solution, u{x,t ) = F(x~Vt) , of the nonlinear Schrodinger 

2 

equation: \U t +(XU xx + j3\ll\ U = 0 (where (X and are real constants). Hence show that there is a solution 
F = a Cxp[i/t(x — fi7)], for a and k real, and determine V. 


x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x- 

********** 
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4 Similarity solutions 


We mentioned at the start of this Notebook that similarity solutions are very special, although without doubt interesting 
and important. The basic idea is introduced by considering the classical example of this type - the heat conduction 
equation - and then we shall make some general observations about the fundamental property that underpins the method. 
Also, as we have done in previous chapters, we will consider both simple, linear, classical partial differential equations and 
some nonlinear examples. Finally, we will outline how these special solutions can be extended and then used to obtain 
more general solutions. 

4.1 Introducing the method 

Let us consider the heat conduction (diffusion) equation in one dimension: 

U t = ku xx ( k > 0, constant), 

and seek a solution u(x, t) = F (xt n ) where n is a constant to be determined; we must also expect that this form of 
solution is valid only for certain functions F. We therefore have 

u t = nxt n ~ l F'(rj) and u xx = t 2n F"(rj ), 

where T] = xt n is the similarity variable; the equation becomes 

nxt n ~^F'{rf) = kt 2n F"(rf). 


Now for this to be a consistent and meaningful equation, it must be an ordinary differential equation for F( /)). First, let 
us introduce 7 7 in place of, for example, x i.e. set X = 'll/1 11 . (We could equally replace all the t- dependence by writing 
t = ( 77 /x) 1 / n ,, but our choice is the simpler route.) Thus we obtain the equation in the form 

mit~ l F'(r/) = kt 2n F"(rj) 


and this can be made consistent only if 2/7 = — 1, so n = —1/2 (and we should note that the equation for F poses serious 
difficulties on t = 0 - which was not evident in the original partial differential equation - so we must use t > 0). With 
this prescription for n, we obtain 


F " = — F ' and so F ' = A exp 

2 k 


where A is an arbitrary constant; finally we have 


77 

f 

77' 2 

Jexp 

4k 


d//. 


However, it is slightly neater to change the integration variable by writing lj = 2 ■4ky to give 


4 k 

/ 

c/Vt 

1 

j exp 

ZL 

4k 
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/ 2Vkt 


u( x, t) = B J exp(-y 2 )dy, 


where 5 is a redefined version of A, and a second arbitrary constant is implied by the indefinite integral. 


To see how special all this is, let us suppose that we wished to use this method to find the solution for the temperature, 
u(x,t ) , in an infinite rod ( X > 0 ), described by 

U t = kll xx with u( 0, t) = lift, 77 —> 77^ as X —^ co , for all t > 0 , and u(x, 0) = 77 2 , 

where 7/n , U\ and 77 2 are constants. First, we note that all the given temperatures are constant, whereas the most general 
problem with this configuration would be 

Uq = M()(0, lt\ = U\(t) and 77 2 = u 2( x ) > 

so we have already simplified the problem. Further, in terms of the similarity variable ( 77 = x/*Jt ), we see that 

Tj = 0 on X = 0 for t > 0 

and 77 —> 00 as X —> 00 for any finite t > 0 , 

but that 77 is undefined on t = 0 . We can proceed, therefore, only by replacing the given initial-boundary-value problem 
by the following: 


77(0, t) = Uq for t > 0; 

7/ —> 77 1 as X —> 00 for any finite t > 0; 

77 —> lh as t —> 0 + for X > 0 . 

With this interpretation, the first condition is unchanged i.e. !j = 0 on X = 0 (/ > 0), but the second and third conditions 
are now equivalent: 

fas x^oo fort>0 
[as t —> 0 + forx>0. 

This in turn requires that, for a similarity solution to exist, we must have U\ = llj , so very special indeed; with all these 
choices, we may formulate the problem for 77 = F{ /)): 
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Tj/2^fk 

u(x,t) = ?{rj) =B J exp(-y 2 )dy, 

with -F(O) = lift and F —> ll\ as TJ —> oo. The solution is therefore 


u( X, t) = U 0 + ( U! — U 0 ) 


^X /2 Vkt ^Xp(_ y 2)dy 

J 0 "exp(-y 2 )cV 


which, with 


j:exp(-y 2 )dy 



, is more neatly written as 


xjl^fVX 

(ui-u 0 ) J exp(-y 2 )dy. 
o 

Before we investigate an important underlying property of partial differential equations that, when satisfied, implies the 
existence of a similarity form, we consider, as examples, our other two elementary, classical equations. 

Example 16 

Seek a similarity solution, in the form u(x, t) = F |xt j, of the wave equation U tt — C^U XX = 0 (where C > 0 is a 
constant). 


u( X, t) = u 0 + 


With U(x,t) — F (xt j (although we can note that a function of tx n would work equally well - indeed, the symmetry 
in x and t might suggest something more), we obtain 

u xx =F n F" and u tt = n(n — X)xt n 2 F'+« 2 x 2 f 2/? 2 F", 

which gives 

n(n -1 )xt n ~ 2 F' + n 2 x 2 t 2n ~ 2 F" - c 2 t 2n F" = 0 . 

/ n 

We introduce X = 1J/ t , and so 

n(n - \)r/t~ 2 F' + n 2 T] 2 t~ 2 F" - c 2 t 2n F" = 0, 
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which is consistent only if n = — 1 (and then t ^ 0 ); thus we have 


f 77 2 -c 2 )f " + 2 / 7 F' = 0 and so [ri 1 -c 2 )f ' = A, 

where A is an arbitrary constant. Hence 


f<„=Aj(_L__L 
2c J { n-c 77 + c 


dr/ = — In 
2c 


71-C 


7]+ C 


+ B„ 


T]-C 7] + ( 

where 5 is a second arbitrary constant. This solution, with 71 = xjt , can be expressed as 

X-Ct 


u(x,t ) = Cln 


x + ct 


+ B = C[ln|x - ct\- ln|x + ct |] + B , 


where C has replaced A/2c ; the essential symmetry in x and t is now evident. Indeed, this solution is no more than a 
special example of d’Alembert’s general solution: 


u = F(x - ct) + G(x + ct) . 
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Example 17 

Seek a similarity solution, in the form u(x, v) = F(ycy ), of Laplaces equation U xx + Uy V = 0 . 


This is virtually identical to Example 16; we therefore obtain 

y ln F” +n(n - \)rjy~ 2 F' +n 2 rj 2 y~ 2 F" = 0 
and then, as before, we require n = — 1: 


F " + 2 / 7 F' + 77 2 F = 0 and so (l + / 7 2 )f ' = A 


which gives F(rf) = A arctan( rj) + B , 

where A and B are arbitrary constants. Thus we have the solution 


u(x, y) = Aarctan 




y J 


+ B . 


Comment: It is clear, by virtue of the symmetry in x and y, that another solution is 


u(x, y) =Carctan 


'-1 
V xj 


+ D 


where C and D are arbitrary constants. That these two solutions are the same solution is clear: 


tan 


u- B 
A 


x 

— and so cot 

y 


u- B 


V 


A 


y 

1 

X 


cot 


( u-D k/ 


f U ~ D 1 


Cy} 

+ 

= tan 


i.e. u = Carctar 

— 

l c 2 ) 


l c J 




+ D 


In other words, a choice of the arbitrary constants in the first ( A = —C ,B/ A = —(F)/C + 7r/2 j) recovers the second. 

4.2 Continuous (Lie) groups 

In order to motivate these ideas, let us return to the heat conduction equation: 

Uf — /v l ^ xx , 
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and transform the variables according to 


u = aU , t = J3T, x = yX, 

where Ot , f:> and y are real, non-zero constants. We now regard U = U ( X , T) , so that the equation for U becomes 

a TT l a TT 

— U T -k^rUxX’ 

P y 

2 

and then for the choice /? = y , and any Cl (^ 0), we recover the same equation but written in the new notation: 


U t — kUxx ■ 


We say that the equation is invariant under the transform 


U = aU , X = yX, t = y T usually expressed as U —> au , X — > yx, t —> y"t, 


where ‘ —> ’ means ‘replace the old by the new’, avoiding the need to change variables (and this reinforces the property of 

leaving the equation invariant). Note that this holds for arbitrary, real, non-zero OC and y ; further, because the equation in 

u is linear, the use of Cl here is redundant, so it is sufficient to state that the equation is invariant under the transformation 
2 

X — )y X , t —> y t . In nonlinear equations, however, CC will normally play a critical role, as we shall now see. 
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We consider the Korteweg-de Vries equation 


U t + UU X + U xxx = 0, 

and transform exactly as above, using CL, f:> and y : 

cl a 1 CL 

-u t +—uu x +-^u xxx = 0. 

P Y y 5 

The equation is therefore invariant under the transformation 

—2 3 

a = y and P = y . 

What we have described are two examples of a continuous or Lie group, which we will now define more precisely. [Marius 
Sophus Lie, 1842-1899, Norwegian mathematician who, with Felix Klein, developed the notion of groups in geometry. He 
also discovered contact transformations, and applied the ideas to integration and the classification of partial differential 
equations.] 

Let the transformations that we have just introduced be represented by 71, (for real y V 0 ); further, we consider the 
application of successive transformations: Ty , Tg, 71., etc., and examine the consequences. 

1. We apply Ty followed by Tg ; the result is the same as applying Tyg , so we immediately have the 
multiplication law: TyTg = Tyg. 

2. Now we apply Tg followed by Ty ; this is also the same as applying Tyg (= Tgy ), so the multiplicative law 
is commutative: TyTg = TgTy . 

3. We can now apply a third transformation, 71., to produce the equivalent of applying Tyg s . 

4. From (3) we see that the multiplication law is associative: 

Te ("TjT y j = TgT yg = T yg e = Tg £ Ty = (T f Tg jTy. 

5. There is an identity transformation,: T\ i.e. 7] 71, = 71, 7] = 71,. 

6. There exists an inverse of Ty , namely, Tyjy , because Ty Tyjy = Ty and also Tyjy Ty = Ty , so Tyjy is both 
the left and right inverse of Ty . 

The conditions and properties (l)-(6) show that the elements of Ty , for all real y 0 , form an infinite, continuous (Lie) 
group, where y is usually called the parameter of the group. 

The property that an equation is invariant under such a transformation suggests that solutions of the equation might be 

constructed that possess the same property i.e. the solution is also invariant under the transformation. For example, the heat 

2 

conduction equation (see above) is invariant under the transformation X —> yx , t —> y t, so we might seek a solution 
that is a function of jt or, equivalently, xjpt - which is exactly what we found in §4.1. Similarly, the wave equation 
and Laplace’s equation are clearly invariant under the transformation X —> yx , t —> yt and X —>■ yx , y —> yy, 
respectively, so solutions that are functions of x/1 and x/y, respectively, are to be expected; see Examples 16 and 17. 
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Example 18 

Find the transformation that leaves ll t +U XXX = 0 (the linearised Korteweg-de Vries equation) invariant, and hence 
seek an appropriate solution (but do not solve the resulting ordinary differential equation). 


Given U t + U xxx = 0 , we transform according to 


t —» at , X —» /?X and then a \l 


+ /? ? \l x 


= 0, 


(and we do not need to do anything more because the equation is linear); so we must choose OC = 
/ 1/3 

more conveniently, x t ' . Hence we seek a solution 


P 


3 i.e. use x /if 


or, 


u(x,t) = F(x/t 1 / 3 ) 


which gives 

-^xr 4/3 F , + (t“ 1/3 ) 3 F ,,, = 0 or F - -1 r/F ' = 0 (7/=x/t 1/3 ). 


Comment: this ordinary differential equation can b e solved; we set F G/) — G (/// 3 j, which gives G ~CG = 0 

/ 1/3 

where C, = 77 / 3 ' ; this is Airy’s equation with a solution G = Aj(<Q (which is related to Bessel functions of order 
1/3 ). Thus we have a solution for u: 


x/(3t) V 3 

u(x, t) = C J Aj(^) 

where C is an arbitrary constant. 


4.3 Similarity solutions of other equations 

In order to demonstrate how this approach can be extended, let us consider, for example, the higher-dimensional heat 
conduction (diffusion) equation 


Ut = k(u xx + u yy ); 


see §2.3 and Example 13. First we transform according to 

t —> at , x —> fix , y —» yy 
(and we note that the equation is linear), and so we obtain 
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a % = 2 u XK + y 2 u yy )„ 


2 2 

which leads to the choice a = fi~ = y . Let us therefore seek a solution using a pair of similarity variables: 


u(x, y,t) = F(x/Vt,y/Vt)„ 

which gives --1 xt _3/ 2 F^ - ^y t -3 / 2 F^ = k(t _1 F w +t _1 F^, 

where 1J = xj yft , ^ = yj yft . This equation can be simplified (for t > 0), so that 

^) =k ( F w +F ^)'> 


which we may solve in any suitable way - by separation of variables, for example! We set F{lj= E( lj) X(^): 

1 

or — 

2 

so we introduce kE " + f >jE ' = AE , where A is the separation constant, leaving kX" + f ^X' = —AX. One simple 
solution is available if we choose A = 0, for then 


-jr/EX -l£EX ' = k(E"X + EX") 


E' 1 -X' 

n -£-: 

E X 


f r " 




X 


"h 


X 


E </?) = Aj'Txpf //' 2 /4k jd//'; X (A) = B exp(-f' 2 /4k)df'. 
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and so we obtain 


u(x, y.fl = c(j^«p(V)d7|j y/2jtF «p(-^)df).. 


where A, B and then C are arbitrary constants. 

Finally, we consider three equations that, even though we have discussed them previously, warrant further investigation 
under this new heading. 

Example 19 

Seek a similarity solution of the Burgers equation U t + Ull x = kll xx (see §3.3), k > 0 constant, that decays as |x| —> oo . 


First we transform U —> CCU , t —> (5t , X —> yx , to give 



Thus we seek a solution u( X, t) = —j= F (x / Vt j, where we have interpreted the choice of (X ( = y ' = fj ^ _ ) in 
terms of t rather than x, because the equation involves only one derivative in t; this yields 



which, for t > 0 , becomes 


2 FF '-(F + 7f ') = kF " ( 77 = x/Vt). 


This equation can be integrated once, to give 


F 2 - r/F = kF' + A , 


where A is an arbitrary constant. However, given that F —> 0 (and we shall assume that F' does likewise) as 1 77 —> co 
, we must have A = 0. Thus we now have 


F 2 - T)F = kF’, 


which is a simple Riccati equation; to solve this, we set F = /J (f)' j(f) , where // is a constant to be chosen, then 
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and we select fl = —k. We are left with 

k 0" = -rj</)' i.e. <p(rf) = Aj^exp(-//' 2 /4k)d7/'; 


see §4.1 (again, A is an arbitrary constant). Thus we have a solution 

k exp(-x 2 /4kt) 


u(x,t) =- 


Vt 


B+\* ,2 ^exp(-y 2 )dy 


where B is an arbitrary constant (obtained by electing to define the lower limit of the integral and then dividing throughout 
by A). 


Our next example also refers back to one of our important earlier equations, discussed in Example 14. 

Example 20 

Find the ordinary differential equation associated with the similarity solution oftheKdV equation U t +6 UU X +U XXX = 0 . 


We transform as in the previous example, to produce 


j3 \i t + 6ay L uu x + y 3 u xxx = 0 and so a = y z , j3 = y 


-3 


-2 


Thus we seek a solution ll( X, t) = 


1 


f \ 

X 


(3t) 


1/3 


(3t) 


1/3 


,, where the ‘3’ is purely a convenience; we obtain 


-^3~ 2 / 3 r 5 / 3 F-(3ty 2 / 3 x(3ty 4 / 3 F' + 6(3ty 5 / 3 FF' + (3ty 5 / 3 F'" = 0 


For t > 0 , this simplifies to give 

F'" + (6 F - rj)F' - IF = 0 where ? 7 = jc/ (30 1/3 , 

which is the required equation. 


Comment: This equation can be transformed into a Painleve equation of the second kind, about which much is known. 
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Example 21 

Find a similarity solution of Uj- = (uUx ) (cf. Example 15) which satisfies u( X, t) dx = Constant and which is 
zero at a finite, non-zero value of the similarity variable. 


The transformation (as used above) gives 

P~ y \\ = ay~ 2 ( uu x ) x and so/? = y 2 /or, 

on the basis of this, the most general form that the similarity solution can take is 

u(x, t) = t n / (2-n) F (x/1 1 ^ 2-1 ^ ), 

for any constant n (^ 2 ). Now let us examine the integral constraint, for any t > 0: 

OO oo 

JuUt)dx = t n/(2-n) J F(x/t 1/(2_n) )dx 

—oo —oo 

oo 

= t (n+i)/(n-2) J p( 77 ) d 77 (where 77= x/t 1 ^" -2 ) ), 
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which is a constant only if (n + \ )/(1 — tl) = 0, so n = — 1, producing the similarity form 

u(x,t)=r 1 / 3 F(x/t 1 / 3 ). 

Thus we obtain 

-±r 4 / 3 F 

or (for f > 0) (FF ')' + f(F +lf ') = 0 
which integrates immediately to give 

3 FF' + rjF = A, 

where A is an arbitrary constant. Now for F = 0 at finite // (assuming that F’ is finite at this same point) requires 
A = 0, so we obtain for every point where F A () ■. 

3F ' = -rj and so F (77 ) = ^(b - 7/ 2 ) 

which is zero at T] = , where B > 0 is another arbitrary constant. (We can see that F ^ VB j is finite.) Thus 

u( x, t) = t-V 3 ( B - x 2 /t 2 / 3 )/6. - 

4.4 More general solutions from similarity solutions 

We describe this important idea by analysing a classical example of this type: the heat conduction (diffusion) equation 

U t = kll xx ( k > 0, constant). 

This equation has a similarity solution (see §4.1) of the form 

u(x,t) = F (x/Vt) ; 

further, the equation is linear, so the transformation U —> CCU plays no obvious role (but note what we did in Example 
21). Let us investigate the consequences of seeking a more general solution, for some F • 

u( X, t) = t n F (x/ Vt) , 

which we should expect is permitted because of the freedom to use any a in the group transformation. Thus we obtain 
nt n ~ l F-^xt n ~^ 2 F’ = kt n ~ l F" or nt n ~ l F-^j]t n ~ l F' = kt n ~ l F" 
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which, for t > 0 and any n, simplifies to 


nF-jjiF' = kF". 

It is convenient to write F(rf) = e m1 ^ G( /]) , for some constant m, and then 

nG ~^ rj (G '~2vnr/G) = k^G"-4m^G / -2mG + 4m 2 ^ 2 G) 

where we choose m = 1/4 k and elect to work with the case n = — 1/2 (and note that we may choose any n, and then 
investigate the solution for G). The equation for G now becomes 

^r/G' = kG" and so G' = Aexp(/7 2 /4k), 

where A is an arbitrary constant. However, this solution implies that F = e 1 G will not decay at infinity - indeed, in 
general, it will grow linearly in ?] ; thus we must choose A = 0 and then G = constant = B . Thus we have the solution 


u(x, t) =-^=exp(-x 2 /4kt). 


This new solution, which is an extension of the classical similarity solution, has an important property. Let us examine 

b 


J u( x, t) dx = -j= Jexp(- x 2 /4kt) dx, 


which corresponds to the calculation that we presented in Example 21; we change the integration variable (X = 2 sfkty ) 
to give 

OO 0° 2 

Ju(> 4 t)dx = 2Bv/k^ Je -y dy = 2B- N /^ : k, 

—OO —OO 

which is a constant i.e. not a function of time. How can we exploit this? 

First, let us observe that the original equation, U t = ku xx , is unchanged if we perform a constant origin shift (in x or 
f); for example, on setting X = a + X , where a is a constant, the equation becomes 

U t = kuxx ■ 

B 


Thus, given that U = exp|— X 2 /4ktj or, rather, — ^=exp|— X 2 /4ktj is a solution, then 


so is 


, x B 
u( x, t) =- r exp 

Vt 


-(x-a)74kt 


But then another solution is 


B r 

u(x,t) =-^=Jexp 


-(x-a)74kt 


da, 
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because the original partial differential equation is linear with constant coefficients, and so we may apply an integral 
operator to the equation, providing the integration is not with respect to either x or t (and the integral exists, of course). 
By the same token, yet another solution is 

1 


u(x, t) = —r= f B(a)exp -(x-a) 2 /4kt da 

Vt J L J 


at least for suitable B(a ) . 


Let us therefore examine the solution (which, if it is thought necessary, can be verified as such by direct substitution into 
the equation): 


1 r 

u(x, t) =-j= J B(a)exp 


(x-a)‘ 

4kt 


da ; 


in particular we ask what value this takes on t = 0 . To accomplish this, it is most natural to change the integration 
variable; for any t > 0, we write a — X = 2~Jkty and^so obtain 

u(x, t) = 2x/k J B(x + 2A/kty) e _y dy. 


This is readily evaluated on t = 0 , and yields 

OO 

u(x,0) = 2x/kB(x) |e _y dy = 2 A /k^rB(x), 


which demonstrates that our new solution, based on a similarity structure, provides the solution to the general initial- 
value problem. That is, given u(x, 0) (with suitable conditions at infinity), the solution is as just described, provided 
that the resulting integral exists (which will very often be the case because of the strong exponential decay at infinity). We 
see that we have evaluated the transformed integral at t — 0 , which may seem a slightly dubious manoeuvre (because 
the change of variable requires t > 0); a more careful analysis adopted for t —> 0 + , using more sophisticated methods, 
confirms the correction of our result. 

Example 22 

Use the result in Example 18, and in the comment following it, to find the solution of U t + U XXX = 0 which satisfies 
«(x,0) = /(x). 


The conventional similarity solution can be expressed as 

u(x, t) = F (x/t 1 / 3 ), 


but because the equation is linear, we may seek a generalised version of this: 

u(x,t) = t n F (x/t 1/3 ). 
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Thus we get 

nt n ~ l F - 1 x t n ~ 4 / 3 F' + r l t n F'" = 0 

which simplifies to nF — - j 1)F' + F"' = 0; 


guided by the introductory example (so on this basis we would hope to find that n = —1/3), we try a solution 
F = A j ( /// 3 1 3 j , where 1j = xjt ' 3 and A j (£) is a solution of the Airy equation: w" — t^W = 0 . Thus we obtain 


nA i -13 _1/ 3 r/k [ + 1A '"= 0 and so A {Qk x )' = (A g\ x )' = 0, 


if we select n = —1/3 (as expected). Hence we have a solution 

u(x -= ^w A i( x /( 3tV3 )) 

where C is an arbitrary constant. 


But we may replace x by X 
solution: 


a (an origin shift), and then write C = C(a ) ; integration over all a then yields another 



f \ 

J C(a)Aj 

—oo 

x-a 

l> V3 )J 


da . 
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1/3 

We now introduce x — a = (3t) y , to give 

OO 

u( x, t) = 3 1,/3 J C (x - (3t) 1//3 y)A j (y) dy 

—00 

and then evaluation on t = 0 yields 

OO 

u(x,0)=3 1 / 3 C(x) J A; (y) dy. 

00 —00 

Let us write J A j ( y) d y = k (it turns out that k = 1), then we require 

—OO 

3 1//3 £C(x) = f(x) 


and so the initial-value problem has the solution 


OO 

u(x,t) =—1 f (a)A 
(3t) 1/3 ki 


f \ 

x- a 

1 - ' 1/3 


.(*) 


da 


vv-'v y 


Exercises 4 

1. Construct the continuous group transformation for the equation U t + IIU x = tu xx , and hence seek an appropriate 
similarity solution (by deriving the ordinary differential equation satisfied by this solution). 

I |2 

2. Repeat Q.l for the nonlinear Schrddinger equation iu t + U xx + \u\ U = 0. 

X-jfX-X-X-X-X-X-jfX-X-X-X-X-X-X- 
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Answers 

Exercises 2 

1 . u(r,t) = R(r)T(t): T" + CQ~c~ T = 0, rR" + R' + a?' R = 0 (a Bessel equation). 

_i 2 

2 . u(r,t ) = r v{r,t ) : v tt — c v xx = 0 - the classical wave equation. 

3. u(x,t) = X(x)T(t): u = (C + Xx)/(l + At) . 

Exercises 3 

1 . u(x,t) = F(x — Vt) : kF" + VF' + AF = 0 with condition V 2 < 4/1 k" . 

2 . u(x,t) = F(x -Vt ): -iVF' + aF" + fllFl 2 F = 0 and then V = ak -P-CL 1 . 

k 

Exercises 4 

1 . u(x,t) = F(x/t): F" +(rj-F)F' = 0 where 77 = x/t . 

-i 

2 . u(x,t) =—. (x/Vt): F "-li(/)F )'+|F | 2 F = 0 where 77 = x/Vt. 

V t 

**************** 

********** 
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nonlinear diffusion 
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parabolic 

104.109 


quasi-linear 

13,15,53 


quasi-linear (n) 

23 


semi-linear 

13,22,54 


shallow water 

88,92 


sine-Gordon 

119 


Tricomi 
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wave 
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wave (solution) 
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125 

Euler's equation 


84.113 
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explicit solution 


61 

First-order 

equation 

55 

first-order, linear 

equation 

56 
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equation 
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functions 

class Cl 

15 

General equation 
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general solution 


36 
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similarity solution 

132 
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125 


with initial data 

146 

hodograph transformation 


85,90 

hyperbolic 

equation 

104.109 

hyperbolic case 


74 


canonical form 

77 

hyperbolic equation 

Cauchy data 

81 


Goursat problem 

81 

Implicit solution 


61 

initial conditions 


81 

initial data 


65 


for heat cond. Eqn 
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general equation 

34 

initial/boundary data 
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82 
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19 

initial-value problem 


18 

invariant 


137 

invariants 

Riemann 

94 

Jacobi, KGJ 


71 

Jacobian 


19,37,71,72,90 

Kink solution 


121 

Korteweg-de Vries (KdV) 

similarity solution 

142 

eqn 

travelling wave 
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Lagrange, J.-L. 


13 

Laplace's equation 
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plane polars 
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separation of variables 
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similarity solution 

136 
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Laplacian 
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Lie 

continuous 

136 

Lie group 


136 

linear equation 


13,22 

linear Korteweg-de Vries 

similarity (extended) 

146 


similarity solution 

139 


with initial data 

147 

Mach lines 


87 

Mach number 


84 

Mach, E 


84 

Monge cone 


26,27,29 

Monge, G. 


13,26 

Neumann problem 


83 

Neumann, CG 


83 

nonlinear diffusion equation 

separation of variables 

118 

Parabolic 

equation 

104.109 

parabolic case 


74 


canonical form 

77 

parabolic equation 

initial/boundary data 

82 

parameter of the group 


138 

plane polars 


112 

Quasi-linear 

equation 

53 

quasi-linear equation 


13,15 

quasi-linear equation n 

variables 


23 

quasi-linear, first order 

equation 

59 

Rectangular membrane 

vibrations 

116 

Riemann invariants 


94 
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94 

Riemann, GFB 


94 
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83 

Robin, G 


83 
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73 

self-similar 
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semi-linear 

2nd order equation 

71 


equation 

54 

semi-linear equation 


13,22 

separation constant 


107 

separation of variables 

cylindrical coords 

116 


cylindrical symmetry 

110 


heat conduction eqn 

109 


in 2D rect. Cart. 

114.115 


in 3D rect. Cart. 

125 


Laplace's equation 

109 


nonlinear diffusion 

118 


plane polars 

112 


sine-Gordon eqn 

119 


spherical coords 

117 


wave equation 

107 

shallow-water equations 


88,92 


Riemann invariants 

94 


simple waves 

96 

shock waves 


87 

shock/discontinuity 


97 

signalling problem 

wave equation 

70 

similarity solution 

Burgers equation 

141 


heat conduction eqn 

132 


heat conduction in 2D 

139 


Korteweg-de Vries eqn 

142 


Laplace's equation 

136 


linear KdV 

139 


wave equation 

134 

similarity solution 

heat conduction eqn 

144 

(extended) 

linear KdV 

146 

similarity variable 
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simple waves 

shallow water eqns 

96 

sine-Gordon equation 

separation of variables 

119 

singular solution 


36,37 

solitary wave 
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sine-Gordon 
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soliton problem 
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solution 

d'Alembert's 
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84 

speed of propagation 


66 

speed of sound 


84 
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Sturm-Liouville problem 


108 

subsonic 


84 

supersonic 


85 

Tangent plane 


15,26 

tangent vector 


26 

Taylor shock profile 


128 

transonic 


85 

travelling wave 

Burgers equation 
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heat conduction eqn 
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heat conduction in 2D 

125 


Korteweg-de Vries eqn 

128 


Laplace's equation 
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of wave equation 
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Tricomi equation 
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85 
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Wave equation 
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Cauchy data 
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domain of dependence 
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68 
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separation of variables 

107 
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similarity solution 
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